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A. 


REVIEW  OF  OCEAN  ENVIRONMENT 


The  ocean  envirorment  Is  characterized  by  waves,  wind  and  current. 
The  waves  are  typically  Irregular  (confused  or  random  seas).  Some 
waves  are  generated  locally  by  the  wind,  and  some  waves  are  generated 
great  distances  away.  The  wind  Is  unsteady,  with  gusts.  The  wind 
varies  with  height  above  water.  The  current  Is  caused  by  the  wind, 
by  waves,  by  the  tide,  and  by  global  temperature  differences.  The 
current  varies  with  depth.  All  of  these  characteristics  vary  with 
time. 

A. I  IRREGULAR  WAVES 


Irregular  waves  (a  random  sea)  can  be  described  as  the  sum  of  an 
Infinite  number  of  Individual  regular  (sinusoidal)  waves  of  different 
amplitude,  frequency,  and  phase  (Figure  A-1).  Therefore,  the 
randomly  varying  sea  surface  elevation  can  be  represented  by  a 
Fourier  series. 


N 

n(t)  =  E  aTcos(a>.*t+*.) 
1=1  ^ 


where  ri(t)  1s  the  water  surface  elevation  measured  from 

still  water  level. 


a^  Is  the  amplitude  of  each  component  regular  wave, 

u>^  Is  the  frequency  of  each  component  regular  wave. 

Is  the  phase  angle  of  each  component  regular 
wave,  and 


t  Is  time. 
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The  most  distinctive  feature  of  a  random  sea  Is  that  It  never  repeats 
Its  pattern  and  it  Is  Impossible  to  predict  Its  shape.  -Therefore, 
total  energy  Is  used  to  define  a  particular  sea.  The  energy  (E)  In 
an  Individual  regular  wave  per  unit  surface  area  Is, 

E  =  i5*p*g*a^ 

and  the  total  energy  of  the  sea  Is  the  sum  of  the  energies  of  the 
constituent  regular  waves. 


E  =  J5*o*g*  I  aj 
1*1  ^ 


The  total  energy  of  the  sea  Is  distributed  according  to  the 
frequencies  of  the  various  wave  components.  The  amount  of  energy  per 
unit  surface  area  within  the  small  frequency  band  (w^,  w^+dw)  Is, 

E(w^)  *  %*p*g*a^^*du 

The  total  energy  of  the  sea  Is  then  the  sum  of  the  energies  within 
the  Individual  wave  components.  If  the  sea  Is  made  up  of  an  Infinite 
number  of  waves,  the  energies  of  the  waves  form  a  smooth  curve,  and 
the  above  summation  may  be  replaced  by  an  Integral. 


E  *  p*g*  J”  J5*a^*du 


The  smooth  distribution  of  the  wave  energy  Is  called  the  energy 
spectrum  or  wave  spectrum  of  the  random  sea,  and  Is  often  designated 
as  S((ii).  A  wave  spectrum  Is  normally  depicted  as  a  curve  with  an 
ordinate  of  energy  and  an  abscissa  of  frequency.  A  typical  wave 
spectrum  has  a  central  peak  with  a  tapered  energy  distribution  either 
side  of  the  peak. 
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The  recommended  form  of  displaying  a  wave  spectrum  Is  with  an 
ordinate  of  Js*a^  and  an  abscissa  of  u,  radial  frequency.  However, 
since  the  engineer  will  encounter  wave  spectra  equations  in  a  number 
of  forms,  using  various  bases  and  units,  the  applicable  conversion 
factors  are  provided  in  the  following  sections. 

Spectrum  Basis 

The  recommended  spectrum  basis  is  half  amplitude  squared  or  energy. 
Often  spectrum  equations  having  a  different  basis  are  encountered. 
Before  any  statistical  calculations  are  performed  with  a  spectrum 
equation,  the  equation  should  be  converted  to  the  recommended  basis. 

For  a  "half  amplitude"  or  "energy"  spectrum,  the  basis  is  one-half 
times  the  amplitude  squared. 

S(w)  * 

S(«)da.  *  E  /  (pg) 


where, 

S  is  the  spectral  ordinate, 

u  is  the  radial  frequency, 

n  is  the  wave  amplitude  of  the  constituent  wave  of 

frequency, 

E  is  the  energy  content  of  the  constituent  wave  of 

frequency,  u. 

For  an  "amplitude"  spectrum,  the  basis  is  amplitude  squared. 

S(u)  = 
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For  a  "height"  spectrum,  the  basis  Is  height  squared. 

S(m)  *  h* 

S(«)  - 

where, 

h  Is  the  height  of  the  constituent  wave  of 

frequency,  &>. 

For  a  "height  double"  spectrum,  the  basis  Is  two  times  the  height 
squared. 

S(«)  -  s*h* 

S(«)  -  16*(is*n*) 

The  basis  of  the  spectrum  must  be  determined  before  the  spectrum  1s 
used  In  an  analysis,  because  the  ordinate  of  one  representation  of 
the  spectrum  may  be  as  much  as  16  times  as  great  as  the  ordinate  of 
another  representation. 

Units 

The  spectrum  equation  may  be  expressed  In  terms  of  radial  frequency, 
circular  frequency,  or  period.  Conversion  between  circular  frequency 
and  radial  frequency  is  accomplished  by  multiplying  by  the 
constant,  2ir. 

u  »  2ii*f 

SM  .  S(f)  /  (2w) 
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where. 


f  Is  the  circular  frequency. 

The  conversion  between  period  and  radial  frequency  Is  more 
complicated. 


ui  -  2»/T 
S(f)  •  t'*s(T) 

S(.)  =  t’*S(T)  /  (2.) 


where. 


T  Is  the  period. 

When  converting  between  period  and  frequency,  the  abscissa  axis  Is 
reversed.  Zero  period  becomes  Infinite  frequency,  and  Infinite 
period  becomes  zero  frequency. 

Wave  spectrum  equations  may  be  used  with  any  length  units  by 
remembering  that  the  spectrum  ordinate  Is  proportional  to  amplitude 
squared  or  height  squared. 

The  mathematical  formulation  for  the  wave  spectrum  equation  will 
often  Include  the  significant  height  squared  or  the  gravitational 
constant  squared,  which  when  entered  In  the  appropriate  units  will 
convert  the  equation  to  the  desired  length  units. 

A. 2  PROBABILITY  CHARACTERISTICS  OF  WAVE  SPECTRA 


The  characteristics  of  ocean  waves  are  determined  by  assuming  that 
the  randomness  of  the  surface  of  the  sea  can  be  described  by  two 
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common  probability  distributions,  the  Gaussian  (or  normal) 
distribution  and  the  Rayleigh  distribution.  These  probability 
distributions  are  used  to  define  the  distribution  of  wave  elevations,  n, 
and  of  wave  heights,  H,  respectively. 

A. 2.1  Characteristic  Frequencies  and  Periods 

For  design  purposes  sea  spectra  equations  are  selected  to  represent 
middle  aged  seas  that  would  exist  some  time  after  a  storm,  yet  which 
are  still  young  enough  to  have  a  good  dispersion  of  wave 

frequencies.  The  primary  assumption  about  the  design  seas  Is  that 
the  wave  elevations  follow  a  Gaussian  or  normal  distribution. 
Samples  of  wave  records  tend  to  support  this  assumption.  In 

conjunction  with  the  Gaussian  distribution  assumption,  the  wave 

elevations  are  assumed  to  have  a  zero  mean.  Digitized  wave  records 
tend  to  have  a  slight  drift  of  the  mean  away  from  zero,  usually 
attributed  to  tide  or  Instrument  drift.  The  Gaussian  distribution 
assumption  Is  equivalent  to  assuming  that  the  phase  angles  of  the 
constituent  waves  within  a  wave  spectrum,  are  uniformly  distributed. 

The  Gaussian  distribution  allows  one  to  calculate  statistical 
parameters  which  are  used  to  describe  the  random  sea.  The  mean 
elevation  of  the  water  surface  Is  the  first  moment  of  the  Gaussian 
probability  density  function.  The  mean-square  is  the  second  moment 
taken  about  zero,  and  the  root -mean-square  Is  the  positive  square 

root  of  the  mean-square.  The  variance  Is  the  second  moment  taken 

about  the  mean  value.  The  standard  deviation  Is  the  positive  square 
root  of  the  variance.  Since  the  wave  elevations  are  assumed  to  have 
a  zero  mean  value,  the  variance  Is  equal  to  the  mean-square,  and  the 
standard  deviation  1s  equal  to  the  root-mean-square.  In  present 
practice,  the  area  under  a  random  wave  energy  spectrum  Is  equated  to 
the  variance. 

In  a  similar  way,  the  characteristic  frequencies  and  periods  of  a 
wave  spectrum  are  defined  in  terms  of  the  shape,  the  area,  and/or  the 
area  moments  of  the  wave  spectrum.  Depending  upon  the 
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particular  wave  spectrum  formula,  these  characteristic  periods  may  or 
may  not  reflect  any  real  period.  The  area  and  area  moments  are 
calculated  as  follows. 

Area; 


m^  *  J*^  (i)^*S(u)du 


Nth  Area  Moment: 


m^  *  J*o 


The  characteristic  frequencies  and  periods  are  defined  as  follows. 
Peak  frequency 

The  peak  frequency  is  the  frequency  at  which  the  spectral  ordinate, 
S(u)  Is  a  maximum. 

Tp!  Peak  period 

The  peak  period  Is  the  period  corresponding  to  the  frequency  at  which 
S(u)  Is  a  maximum. 


T„;  Modal  period 
m 

The  modal  period  Is  the  period  at  which  S(T)  Is  a  maximum.  Since  the 
spectrum  equations  In  terms  of  frequency  and  In  terms  of  period 
differ  by  the  period  squared  factor,  the  modal  period  Is  shifted  away 
from  the  peak  period. 


T^:  Visually  Observed  Period,  or  Mean  Period,  or  Apparent 

Period 

The  visually  observed  period  Is  the  centroid  of  the  S(u)  spectrum. 
The  International  Ship  Structures  Congress  (ISSC)  and  some 
environmental  reporting  agencies  have  adopted  T^  as  the  period 
visually  estimated  by  observers. 

T,  ■  2.*(ino/iiij)'« 

T^:  Average  Zero-upcrossing  period  or  Average  Period 

The  average  zero-upcrossing  period  1s  the  average  period  between 
successive  zero  up-crossings.  The  average  period  may  be  obtained 
from  a  wave  record  with  reasonable  accuracy. 

T^  :  Crest  Period 

The  crest  period  Is  the  average  period  between  successive  crests. 
The  crest  period  may  be  taken  from  a  wave  record,  but  its  accuracy  Is 
dependent  upon  the  resolution  of  the  wave  measurement  and  recording 
equipment  and  the  sampling  rate. 

T^;  Significant  Period 

The  significant  period  Is  the  average  period  of  the  highest  one-third 
of  the  waves.  Some  environmental  reporting  agencies  give  the  sea 
characteristics  using  T^  and  H^,  the  significant  wave  height.  There 
are  two  equations  relating  T^  to  Tp. 

Tj  -  0.8568*Tp,  Old 


Tg  -  0.9457*Tp,  New 

The  first  equation  applies  to  original  Bretschnelder  wave  spectrun, 
and  the  second  Is  the  result  of  recent  wave  studies  (See  Reference 
A.l). 

The  peak  period,  Tp,  Is  an  unambiguous  property  of  all  common  wave 
spectra,  and  Is  therefore  the  preferred  period  to  use  In  describing  a 
random  sea. 

A. 2. 2  Characteristic  Wave  Heights 

From  the  assumption  that  the  wave  elevations  tend  to  follow  a 
Gaussian  distribution.  It  Is  possible  to  show  that  the  wave  heights 
follow  a  Rayleigh  distribution.  Since  wave  heights  are  measured  from 
a  through  to  succeeding  crest,  wave  heights  are  always  positive  which 
agrees  with  the  non- zero  property  of  the  Rayleigh  probability 
density.  From  the  associated  property  that  the  wave  heights  follow  a 
Rayleigh  distribution,  the  expected  wave  height,  the  significant  wave 
height,  and  extreme  wave  heights  may  be  calculated.  The  equation  for 
the  average  height  of  the  one-over-nth  of  the  highest  waves  Is  as 
follows. 

H^/n/  (iflo)  **  =  2*  (2*10(0))**  + 

n*(2ir)***{l-erf((ln(n)>**]} 

where: 

m  Is  the  variance  or  the  area  under  the  energy 

0 

spectrum. 

In  Is  the  natural  logrithm, 

erf  Is  the  error  function,  (the  error  function  Is 

explained  and  tables  of  error  function  values  are 
available  In  mathematics  table  books.) 
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The  characteristic  wave  heights  of  a  spectrum  are  related  to  the 
total  energy  In  the  spectrum.  The  energy  Is  proportional  to  the  area 
under  the  )s*a^  spectrum. 

Average  Wave  Height 

The  average  or  mean  height  of  all  of  the  waves  Is  found  by  setting 
n-1. 

»  2.51*(m^)’* 

Significant  Height 

The  significant  height  Is  the  average  height  of  the  highest  one-third 
of  all  the  waves,  often  denoted  as 

H^  »  4.00*(mJ** 

H,„^:  Maximum  Height 

The  maximum  height  Is  the  1  rgest  wave  height  expected  among  a  large 
number  of  waves,  (n  on  the  order  of  1000),  or  over  a  long  sampling 
period,  (t  on  the  order  of  hours). 

The  maximum  wave  height  Is  often  taken  to  be  the  average  of  the 
1/ 1000th  highest  waves. 

H  ,  -  7.94*(m  )**  -  1.985*H^ 

i/iooo  'o'  S 

Using  the  one-over-nth  equation  and  neglecting  the  second  term  gives 
the  following  equation. 

•  2*1  ln(n)|'» 

or 
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"i/n  *  I'^s 


For  n  «  1000,  this  gives. 


H 


l/lOOO 


For  a  given  observation  tine,  t.  In  hours,  the  nost  probable  extreae 
wave  height  Is  given  by  the  following  equation. 

“max  •  2*(2*»^*ln(3600*t/T2)l^ 


The  3600*t/T2  is  the  average  nuaber  of  zero  up-crossings  In  time,  t. 
A.3  WAVE  SPECTRA  FORMULAS 


The  Bretschnelder  and  Plerson-Hoskowitz  spectra  are  the  best  known  of 
the  one-dlaenslonal  frequency  spectra  that  have  been  used  to  describe 
ocean  waves.  The  JONSWAP  spectrum  Is  a  recent  extension  of  the 
Bretschnelder  spectrum  and  has  an  additional  term  which  may  be  used 
to  give  a  spectrum  with  a  sharper  peak. 

A.3.1  Bretschnelder  and  ISSC  Spectrum 

The  Bretschnelder  (Reference  A. 2)  spectrum  and  the  spectrum  proposed 
as  a  modified  Plerson-Moskowitz  spectrum  by  the  Second  International 
Ship  Structures  Congress  (Reference  A. 4)  are  Identical.  The 
Bretschnelder  equation  In  terms  of  radial  frequency  Is  as  follows. 

S(«)  -  (5/16)*(H5)**(V/«‘)*  exp  (-1.25*(«/«J''‘  I 


where: 


A-ll 


Is  the  significant  wave  height,  and 


(ii  Is  the  frequency  of  maxlniuin  spectral  energy, 

fn 

The  Bretschnelder  equation  may  be  written  In  terms  of  the  peak  period 
Instead  of  the  peak  frequency,  by  substituting  2ir/Tp. 


S(«)  -  (5/16)*(H,)’»  l(^•)‘/(«'*(TJ'•)l* 
expl-l?25*(2.*/(«*T  )*|  >’ 


A. 3. 2  Plerson-Moskowitz  Spectrum 

The  Plerson-Moskowitz  (Reference  A. 4)  spectrum  was  created  to  fit 
North  Atlantic  weather  data.  The  P-M  spectrum  Is  the  same  as  the 
Bretschnelder  spectrum,  but  with  the  and  dependence  merged 
Into  a  single  parameter.  The  frequency  used  In  the  exponential  has 
also  been  made  a  function  of  reported  wind  speed.  The  equation  for 
the  Plerson-Moskowitz  spectrum  Is  as  follows. 

S(u)  »  a*g^/«®)*exp(-6*(u^/u))**  1 


where: 


01  :  0.0081 

8  -  0.74 

*  g/u 

and,  U  Is  the  wind  speed  reported  by  the  weather  ships. 

The  Plerson-Moskowitz  spectrum  equation  may  be  obtained  from  the 
Bretschnelder  equation  by  using  one  of  the  following  relations 


between  H.  and  u.. 

^  m 
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-  0.1610*g/(«„)’ 


or 


»  0. 40125*9/ (H^)^ 

An  Interesting  point  that  may  be  noted  Is  that  If  b  were  set  equal 
to  0.75  Instead  of  0.74,  the  would  be  the  frequency 

corresponding  to  the  modal  period,  T^. 

A. 3. 3  JONSWAP  and  Related  Spectra 

The  JONSWAP  wave  spectrum  equation  resulted  from  the  Joint  North  Sea 
Wave  Project  (Reference  A. 5).  The  JONSWAP  equation  Is  the  original 
Bretschnelder  wave  spectrum  equation  with  an  extra  term  added.  The 
extra  term  may  be  used  to  produce  a  sharply  peaked  spectrum  with  more 
energy  near  the  peak  frequency.  The  JONSWAP  spectrum  can  be  used  to 
represent  the  Bretschnelder  wave  spectrum,  the  original  Plerson- 
Moskowltz  wave  spectrum,  and  the  ISSC  modified  P-M  spectrum.  The 
full  JONSWAP  equation  Is  as  follows. 

S(«)  =  (cij*g^«®)  *  expl-1.25*u./tt|j|'‘*)*T® 


where: 


a  *  exp  1.1s*(m-uiJ*  /  (o*«J^l 

Is  the  frequency  of  maximum  spectral  energy. 

The  Joint  North  Sea  Wave  Project  recommended  the  following  mean 
values  to  represent  the  North  Sea  wave  spectra. 

Y  »  3.3 

a  *  0.07,  for 
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a  =  0.09«  for 

m 

The  value  of  a  1s  found  by  Integrating  the  spectrum  and  adjusting 
w  to  give  the  desired  area. 

The  Bretschnelder  equation  and  the  ISSC  equation  can  be  obtained  by 
setting  the  following  parameter  values. 

Y  »  1.0 

.  .  (5/16)*(Hj.)  =*(„.„)  Vq’ 

The  P1erson>Moskow1tz  equation  Is  obtained  from  the  further 
restriction  that  H.  and  u.  are  related. 

5  ffl 

Hj  =  0.1610*g/(«^)^ 


or 


=  0.140125*(g/H^)’* 


or 


a  »  0.0081 

When  Y  Is  set  to  one  the  JONSWAP  term  Is  effectively  turned  off. 
Without  the  JONSWAP  term,  the  wave  spectrum  equation  can  be 
mathematically  Integrated  to  give  the  following  relationships  among 
the  characteristic  wave  periods. 

Tp  =  1.1362  *  T„ 

Tp  -  1.2957  *  \ 

Tp  «  1.4077  *  T^ 
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Tp  -  1.1671  *  Tj 


For  Y  «  1,  the  fourth  area  moment  Is  Infinite.  The  crest  period, 
Tj.,  Is  therefore  zero. 

For  values  of  y  other  than  one,  the  JONSUAP  equation  cannot  be 
mathematically  Integrated.  The  period  relationships  as  a  function 
of  Ycan  be  calculated  by  numerical  Integration  of  the  wave  spectrum 
equation  over  the  range  from  three- tenths  of  the  peak  frequency  to 
ten  times  the  peak  frequency. 

The  shape  of  the  JONSWAP  spectrum  can  be  further  adjusted  by  changing 
the  values  of  e.  The  e  values  are  sometimes  varied  when  the  JONSUAP 
spectrum  Is  used  to  fit  measured  wave  spectra. 

A. 3. 4  Scott  and  Scott-Wlegel  Spectra 

The  Scott  (Reference  A. 6)  spectrum  was  also  formulated  to  fit  North 
Atlantic  weather  data.  The  Scott  spectrum  Is  the  Derbyshire 
(Reference  A. 7)  spectrum  with  slight  modifications  to  the  constants 
In  the  equation.  The  spectrum  equation  Is  as  follows. 


S(«)  »  0.214*(H5)^*exp(.(«-«J/  {0.065*(u-«^+ 0.26)}**! 


for  -0.26  <  u-«_  <  1.65 
m 


»  0,  elsewhere. 


where 
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Is  the  significant  height. 


=  3.15*r‘-^8,98*r^ 
in 

T  Is  the  characteristic  period  of  the  waves. 

The  U|||1s  the  frequency  of  the  peak  spectral  energy,  but 
unfortunately,  the  period,  T,  used  In  the  equation  for  undoes  not 
correspond  to  any  of  the  mathematical  characteristics  of  the 
spectrum.  The  equation  for  «),gWas  derived  as  a  curve  fit  to  real 
data. 

The  Scott-Wlegel  spectrum  Is  a  Scott  spectrum  modification  that  was 
proposed  by  Wiegel  (Reference  A. 8).  The  constants  are  adjusted  to 
match  the  equation  to  a  ''100~year  storm*  wave  condition.  The  new 
equation  Is  as  follows. 

S(«)  -  0.300*(H5)*exp(.(u-o)^)./ 

{0.0.353*(a.-«^+  0.26)  }  1 


The  Uj^ln  this  equation  Is  1.125  times  that  specified  for  the  Scott 
equation. 

A.4  SELECTING  A  WAVE  SPECTRUM 


Information  about  the  random  sea  characteristics  In  a  particular  area 
Is  derived  by  either  'wave  hindcasting'  or  by  direct  wave 
measurement.  For  many  areas  of  the  world's  oceans,  the  only  data 
available  Is  measured  wind  speeds  and  visually  estimated  wave 
heights.  Sometimes  the  estimated  wave  heights  are  supplemented  by 
estimated  wave  periods.  For  a  few  areas  of  Intense  oil  development, 
such  as  the  North  Sea,  direct  wave  measurement  projects  have  produced 
detailed  wave  spectra  Information. 
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A. 4.1  Wave  Hindcastinq 


Wave  h1ndcast1ng  Is  a  term  used  to  describe  the  process  of  estimating 
the  random  sea  characteristics  of  an  area  based  upon  meteorological 
or  wind  data.  Various  researchers  (References  A. 2,  A. 4,  A. 6,  A. 7  and 
A. 8)  have  attempted  to  derive  a  relationship  between  the  wind  speed 
over  a  recent  period  of  time  and  the  spectrum  of  the  random  sea 
generated  by  the  particular  wind.  The  wind  speed  data  Is  usually 
qualified  by  two  additional  parameters,  the  duration  that  the  wind 
has  been  blowing  at  that  speed  and  the  fetch  or  distance  over  open 
ocean  that  the  wind  has  been  blowing. 

A  set  of  equations  as  derived  by  Bretschnelder  (Reference  A. 2),  which 
relate  wind  speed,  duration  and  fetch  are  as  follows. 

g*H^/U2  =  0.283nanh(0..125*(g*F/U2)°*‘*M 

9*1^  /  (2ir*U)  =  1.2nanh(0.077*(gF/U2)°*'‘* 

g*tn.<n/U  *  6.5882*exp{(0.161*A^-0.3692*A+2.024)’‘ 

+  0.8798*a} 


where 

U  Is  the  wind  speed, 

F  Is  the  fetch, 

A  »  ln(g*F/U^), 

t^^P  Is  the  minimum  duration  for  which  the  fetch  will 

determine  the  significant  height  and  period,  and 

tanh  Is  the  hyperbolic  tangent. 


A-17 


If  the  wind  duration  Is  less  than  tmln,  then  the  third  equation  Is 
used  to  find  the  fetch  which  would  correspond  to  t^^^  »  t. 

For  a  fully  arisen  sea,  the  above  equations  simplify  to  the 
following. 

g*H^/U^  »  0.283 

g*T^/(2ir*U)  =  1.2 

Other  relationships  have  been  developed  1n  the  references.  Often 
specialized  weather/wave  research  companies  have  developed  elaborate 
wave  hindcasting  models  to  derive  the  wave  spectra  characteristics 
for  particular  areas.  However,  the  assumptions  Incorporated  Into 
these  models  have  very  profound  Impact  on  the  outcome. 

A. 4. 2  Direct  Wave  Measurements 


By  Installing  a  wave  probe  or  a  wave  buoy  In  the  ocean  area  of 
Interest,  wave  elevation  histories  may  be  directly  measured.  The 
elevation  of  the  sea  at  a  particular  point  Is  either  recorded  by 
analog  means  or  Is  sampled  at  short  time  Intervals  (typically  one 
second)  and  recorded  digitally.  The  wave  elevations  are  usually 
recorded  Intermittently,  1e.  the  recorder  Is  turned  on  for  say  30  min 
every  four  hours. 

The  wave  records  are  then  reduced  by  computer,  and  the  wave 
characteristics  are  summarized  In  various  ways.  Two  common  ways  of 
summarizing  the  data  are  as  a  wave  scatter  diagram  and/or  as  a  wave 
height  exceedance  diagram. 

The  wave  scatter  diagram  1s  a  grid  with  each  cell  containing  the 
number  or  occurances  of  a  particular  significant  wave  height  range 
and  wave  period  range.  The  wave  period  range  may  be  defined  In  terms 
of  either  peak  period  or  zero-upcrossing  period. 
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The  wave  height  exceedance  diagram  Is  a  curve  showing  the  percentage 
of  the  wave  records  for  which  the  significant  wave  height  was  greater 
that  the  particular  height. 

WAVE  SCAHER  DIAGRAM 


Wave  scatter  diagrams  show  the  occurances  of  combinations  of 
significant  wave  height  and  average  zero-upcrossing  period  over  an 
extended  time  period  such  as  many  years. 

Wave  height  distribution  over  time  can  be  obtained  by  actual  wave 
measurements.  The  heights  and  periods  of  all  waves  1n  a  given 
direction  are  observed  for  short  periods  of  time  at  regular 
intervals.  A  short  time  Interval  of  several  hours  may  be  considered 
constant.  For  this  sea  state,  defined  as  “stationary",  the  mean  zero 
up-crossing  period,  T2,  and  the  significant  wave  height,  H^,  are 
calculated.  The  and  T^  pairs  are  ordered  and  their  probabilities 
of  occurance  written  In  a  matrix  form,  called  a  wave  scatter  diagram. 

Sometimes  wave  scatter  diagrams  are  available  for  the  sea  and  for  the 
swell.  The  sea  scatter  diagram  includes  the  sea  spectra  generated 
locally.  The  swell  scatter  diagram  contains  the  swell  spectra  (or 
regular  waves)  generated  far  from  the  area,  days  before.  Due  to 
greater  energy  losses  In  high  frequency  waves  and  the  continual  phase 
shifting  caused  by  viscosity,  the  energy  in  irregular  seas  tends  to 
shift  toward  longer  periods,  and  the  spectra  becomes  more  peaked  as 
time  passes.  The  energy  in  the  swell  is  concentrated  about  a  single 
long  period/low  frequency,  and  often  the  swell  is  treated  as  a  single 
regular  wave. 

A  typical  wave  scatter  diagram,  presenting  statistical  data  on  the 
occurance  of  significant  wave  height  and  zero  up-crossing  period  per 
wave  direction  is  shown  on  Figure  A-2. 
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Sample  Wave  Scatter  Diagram 


S 

1 

g 

n 

i 

f 

i 

c 

a 

n 

t 

U 

a 

V 

e 

H 

e 

i 

g 

h 

t 

(m) 


12  + - + - + - + - + - - + - + - + - + - + - + 

I  I  I  I  I  I  I  I  i  I  I  i 

11  + - + - + - + - + - + - + - + - + - + - + - + 

I  I  I  I  I  I  I  0.5  I  1.0  I  i  I  I 

10  + - + - + - -* - + - + - + - + - + - - + - + 

I  I  i  I  I  I  I  1.0  I  2.0  I  1.5  I  I  I 

9  + - + - .< - ^ - ♦ - + - + - + - + - ■». - - + 

I  I  I  I  I  0.5  I  1.5  I  2.5  I  3.0  I  0.5  I  I  | 

8  + - + - + - + - - + - + - + - - + - + - + 

I  I  I  I  I  1.0  I  5.0  I  5.5  I  2.5  I  0.5  I  I  I 

7  + - + - + - + - + - + - + - + - + - + - + - -f 

I  I  I  I  I  5.0  113.0  111.0  I  2.0  I  I  I  I 

6  + - - - + - + - + - + - - + - + - - 

I  I  I  0.5  I  6.0  118.0  123.0  I  8.5  I  1.0  |  |  |  | 

5  + - + - ^ - ^ - + - + - + - -f - + - + - + - + 

I  I  I  4.0  126.5  148.5  126.5  I  7.0  |  2.5  I  0.5  I  0.5  I  I 

4  + - 4- — — ♦ - - + - + - + - -f - + - -f - + - + 

I  I  1.5  139.5  179.5  163.5  120.0  I  6.0  I  3.0  I  1.5  I  0.5  I  0.5  I 


0.5 

— - 

150.0 

.+ - + - 

1105.0195.5 

.+ - 

135.0 

.+ - 

111.5 

1  5.5 

.+ - 

1  2.0 

- 

1  1.5 

1 

.•f— 

1 

— + 

1 

--4 

1.5 

159.5 

-+ 

189.0  134.5 

112.0 

1  7.0 

1  4.0 

1  1.5 

1  0.5 

1 

1 

1 

•  •4 

2.5 

118.0 

1  8.0 

1  2.5 

1  2.5 

1  1.5 

.4. 

1  0.5 

1 

1 

.4 

1 

»4.— 

1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

Zero  Up-crossing  Period,  Tz  (sec) 
Sum  of  Occurences  999.5 


Figure  A-2  A  Typical  Wave  Scatter  Diagram  for  the  Central  North  Sea 
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Using  the  significant  wave  height  and  zero  up-crossing  period  from 
the  wave  scatter  diagram  and  selecting  a  representative  sea  spectrum 
formulation,  the  energy  of  each  sea  state  can  be  reconstructed. 

A. 6  WAVE  EXCEEDANCE  CURVE 


A  wave  exceedance  curve  shows  the  number  (percentage)  of  waves  that 
are  greater  than  a  given  wave  height  for  consistent  wave  height 
intervals.  Table  A-1  shows  the  type  of  data  contained  on  a  wave 
exceedance  curve. 


Height  (ft) 

Number  of  Waves 

0 

35,351,396 

5 

3,723,300 

10 

393,887 

15 

41,874 

20 

4,471 

25 

480 

30 

51 

35 

5 

40 

1 

Table  A-1  Wave  Exceedance  Data  for  Campos  Basin 
(Number  of  Waves  from  Northeast) 

This  data  can  be  plotted  on  semi-log  paper  and  closely  approximated 
by  a  straight  line  plot.  Typically,  a  wave  exceedance  H-N  curve  can 
be  defined  with  the  following  equation. 

^  =  ^^m  +  ""z  *  ^09  % 

where 

is  the  maximum  wave  height  for  the  design  life, 
m^  is  the  slope  of  the  H-log  N  curve,  -H^/log  Nj,, 


MIC:500400cc-A 
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is  the  total  number  of  waves  In  the  design  life,  and 

Is  the  number  of  occurances  of  waves  with  height  exceeding 
H. 

A. 7  WAVE  HISTOGRAM  AND  THE  RAYLEIGH  DISTRIBUTION 


Actual  wave  height  measurements  can  be  plotted  to  show  the  number  of 
waves  of  a  given  height  at  equal  wave  height  intervals.  The 
histogram  obtained  can  be  defined  by  a  simple  curve. 

A  simple  curve  that  fits  most  wave  histograms  Is  the  Rayleigh 
distribution.  Past  work  have  shown  that  the  Rayleigh  distribution 
often  allows  accurate  description  of  observed  wave  height 
distributions  over  a  short  term. 

The  Rayleigh  distribution  Is  typically  given  as, 

P(H^)  -  2  *  H^  *  EXP(-H^Vh  *)  *  (1/H 

where 

P(H^)  Is  the  wave  height  percentage  of  occurances, 

H^  Is  the  wave  heights  at  constant  Increments, 

H  ^  1s  the  average  of  all  wave  heights  squared. 

A.8  EXTREME  VALUES  AND  THE  WEIBULL  DISTRIBUTIDN 


For  design  purposes  an  estimate  of  the  maximum  wave  height  (extreme 
value)  Is  required.  The  Rayleigh  distribution  provides  such  an 
estimate  over  a  short  duration.  However,  In  order  to  estimate  the 
extreme  wave  that  may  occur  In  say  IDD  years,  the  Welbull 
distribution  Is  often  used. 
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The  equation  for  the  Welbull  distribution  Is  as  follows. 

P(H)  =  1  -  EXP(  -(  (H.e)/8  )“  1 


where 

P(H)  Is  the  cumulative  probability, 

H  Is  the  extreme  height, 

e  Is  the  location  parameter  that  locates  one  end  of 

the  density  function, 

e  Is  the  scale  parameter,  and 

a  Is  the  shape  parameter. 

By  plotting  the  wave  exceedance  data  on  Welbull  graph  paper,  the 
distribution  can  be  fit  with  a  straight  11ne  and  the  extreme  value 
for  any  cumulative  probability  can  be  found  by  extrapolation. 

A.9  WIND  ENVIRONMENT 

The  wind  environment,  source  of  most  ocean  waves.  Is  random  In 
nature.  The  wind  speed,  its  profile  and  its  directionality  are 
therefore  best  described  by  probabilistic  methods. 


A.9.1  Air  Turbulence,  Surface  Roughness  and  Wind  Profile 


Air  turbulence  and  wind  speed  characteristics  are  primarily 
Influenced  by  the  stability  of  the  a1r  layer  and  terrain.  For 
extreme  wind  gusts  the  influence  of  stability  is  small,  making 
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turbulence  largely  a  function  of  terrain  roughness.  In  an  ocean 
environment,  the  wave  profile  makes  prediction  of  wind 
characteristics  more  difficult.  As  the  wind  speed  Increases,  the 
wave  height  also  Increases,  thereby  Increasing  the  surface 
roughness.  A  surface  roughness  parameter  is  used  as  a  measure  of  the 
retarding  effect  of  water  surface  on  the  wind  speed. 

A  simple  relationship  developed  by  Charnock  (Reference  A. 9)  1s  often 
used  to  define  the  surface  roughness  parameter  and  the  frictional 
velocity  In  terms  of  mean  wind  speed.  Further  discussion  on  surface 
roughness  parameter  and  drag  factor  Is  presented  in  an  ESOU  document 
(Reference  A. 10). 

Full  scale  experiments  carried  out  by  Bell  and  Shears  (Reference 
A.  11)  may  Indicate  that  although  turbulence  will  decay  with  the 
distance  above  sea  surface.  It  may  be  reasonably  constant  to  heights 
that  are  applicable  for  offshore  structures. 

Considering  that  wind  flow  characteristics  are  primarily  Influenced 
by  energy  loss  due  to  surface  friction,  the  mean  wind  profile  for  an 

ocean  environment  may  be  assumed  to  be  similar  to  that  on  land  and  to 

follow  this  power  law: 

*  Vl  “ 

where: 

V  ^2  «  mean  wind  velocity  at  height  z  above  LAT 

V  ,„2i  *  wind  velocity  of  reference  height  above 

UT 

z  >  height  at  point  under  consideration  above 

LAT 
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reference  height,  30  ft  (10  M),  above  LAT 
(typical) 


a  »  height  exponent,  typically  0.13  to  0.15. 

A. 9. 2  Applied.  Mean  and  Cyclic  Velocities 

The  random  wind  velocity  at  height  z  can  be  thought  of  as  a 
combination  of  time-averaged  mean  velocity,  V^,  and  a  time  varying 
cyclic  component,  V2  (t). 

^2  (t)  ■  +  Vj  (t) 

A  range  of  mean  and  associated  cyclic  wind  speeds  can  be  extracted 
from  an  anemogram  and  divided  Into  one-  to  four-hour  groups  over 
which  the  cyclic  component  of  the  wind  speed  Is  approximately 
equal.  By  describing  cyclic  wind  speeds  associated  with  an  average 
value  of  the  mean  component  of  the  wind  speed  over  a  particular 
period  of  time,  a  number  of  pairs  of  mean  and  associated  cyclic 
speeds  can  be  obtained.  In  addition  to  the  applied,  mean  and  cyclic 
wind  speeds  shown  on  Table  A-2,  their  probability  of  occurrence  Is 
necessary  to  generate  a  scatter  diagram.  If  sufficient  data  are  not 
available,  the  number  of  occurrences  can  be  extrapolated  based  on 
similar  data.  Table  A-2  Is  given  only  to  Illustrate  the  wind  make-up 
and  the  uncertainties  associated  with  wind  data. 
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Applied  Wind  Mean  Wind  Cyclic  Wind  Probability 

Speed  Vz(t)  Speed  Vmz  Speed  vz(t)  of 

ft/s  (m/s)  ft/s  (m/s)  ft/s  (m/s)  Occurrence  % 


4.26  (13)  29.5  (9)  13.1  (4)  16.7 
78.7  (24)  62.3  (19)  16.4  (5)  45.8 
101.7  (31)  78.7  (24)  23.0  (7)  12.5 
134.5  (41)  101.7  (31)  32.8  (10)  16.7 
154.0  (50)  124.6  (38)  39.4  (12)  4.2 
180.4  (55)  131.2  (40)  49.2  (15)  4.2 

Table  A-2  Applied,  Mean  and  Cyclic  Wind  Speed  Distribution  for  an 
Extreme  Gust  Environment 


A. 9. 3  Gust  Spectra 


The  power  spectral  density  function  provides  Information  on  the 
energy  content  of  fluctuating  wind  flow  at  each  frequency 
component.  A  study  of  90  strong  winds  over  terrains  of  different 
roughness  In  the  United  States,  Canada,  Great  Britain,  and  Australia 
at  heights  ranging  from  25  feet  (8m)  to  500  feet  (150m)  allowed 
Davenport  (Reference  A. 12)  to  propose  a  power  density  spectrum  of 
along-wind  gust  (the  longitudinal  component  of  gust  velocity). 


A  modified  version  of  the  Davenport  spectrum,  due  to  Harris 
(Reference  ??),  is  given  by: 


nS(n)  ^  4f 

k  V  ^  (2  +  f^)®^* 

m 

where: 
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n 

S(n) 

k 


V 

f 

L 


fluctuating  frequency  2 

power  density  [(m/sec  )/H2l 

surface  roughness  drag  factor  corresponding 
to  the  mean  velocity  at  30  ft  (10m)  (I.e. 
0.0015) 

mean  hourly  wind  speed  at  30  ft  (10m)  m 

non-dimensional  frequency  (nL/V  )  m 

length  scale  of  turbulence  (  1200  to  1800m, 
typical) 


The  Harris  spectra  may  be  used  to  develop  the  wind  spectra  for  each 
one  of  the  mean  wind  speeds  associated  with  the  scatter  diagram. 
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IRREGULAR  WAVES 


Figur*  A-l 


Regular  and  Irragular  Wavaa 
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REVIEW  OF  LINEAR  SYSTEM  RESPONSE  TO  RANDOM  EXCITATION 


B.l  GENERAL 
B.1.1  Introduction 

Spectral  analysis  Is  used  to  determine  the  response  of  linear  systems 
to  random  excitation.  In  the  case  of  offshore  structures,  the  random 
excitation  comes  from  either  Irregular  waves  or  winds.  Typical 
offshore  systems  subjected  to  spectral  analysis  Include  ships, 
semisubmersibles.  Jack-ups,  tension-leg  platforms  and  bottom- 
supported  fixed  platforms.  Responses  of  Interest  Include  motions, 
accelerations,  and  member  Internal  forces,  moments,  and  stresses. 

Floating  units  are  evaluated  by  spectral  analysis  for  motions  In 
random  seas.  The  strength  and  structural  fatigue  Integrity  are  often 
assessed  with  spectral  analysis. 

B.l. 2  Abstract 


A  spectral  analysis  combines  a  set  of  regular  wave  response  amplitude 
operators,  RAOs,  with  a  sea  spectrum  to  produce  a  response 
spectrum.  Characteristics  of  the  response  may  be  calculated  from  the 
response  spectrum,  and  a  random  sea  transfer  function  can  be  derived. 

For  certain  spectral  analyses,  the  sea  spectrum  must  be  modified  to 
produce  a  wave  slope  spectrum  or  to  adjust  the  sea  spectrum  for 
vessel  speed.  A  spreading  function  can  be  applied  to  the  sea 
spectrum  to  model  a  short-crested  random  sea. 

A  wave  force  spectrum  can  be  created  directly  from  force  RAOs  and  the 
sea  spectrum. 

A  regular  wave  transfer  function  Is  found  as  the  solution  to  the 
equations  of  motion.  The  regular  wave  transfer  function  can  be 
expressed  In  terms  of  RAO  and  a  phase  angle. 
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A  white  noise  function  may  be  used  to  represent  a  very  broad  banded 
Input  spectrum.  If  the  response  spectrum  Is  narrow  banded. 

The  extreme  response  can  be  calculated  from  a  given  extreme  wave,  or 
the  extreme  response  may  be  statistically  derived  from  a  set  of 
spectral  analyses. 

The  sea  spectra  used  In  the  computation  of  random  sea  response  can  be 
reduced  In  number  by  selecting  a  smaller  family  of  representative 
spectra,  or  by  creating  a  set  of  mean  spectra. 

B.1.3  Purpose 

The  purpose  of  this  appendix  Is  to  provide  a  background  of  the 
spectral  analysis  method  and  to  clarify  the  concept  of  a  response 
spectrum  and  how  Its  properties  are  derived. 

B.2  RESPONSE  TO  RANDOM  WAVES 


The  spectral  analysis  method  1s  a  means  of  taking  the  known  response 
of  an  offshore  structure  to  regular  waves  and  determining  the 
structure's  response  to  a  random  sea.  The  input  to  the  spectral 
analysis  method  Is  the  response  amplitude  per  unit  wave  amplitude  (or 
equally,  the  response  double  amplitude  per  unit  wave  height)  for  a 
range  of  periods  or  frequencies  of  regular  waves.  These  ratios  of 
response  amplitude  to  wave  amplitude  are  known  as  “Response  Amplitude 
Operators"  or  just  "RAOs."  The  response  of  the  offshore  structure  Is 
first  obtained  for  a  set  of  unit  amplitude,  regular,  sinusicdal 
waves.  The  regular  wave  response  may  be  obtained  either  from  model 
tests  or  from  empirical  or  theoretical  analyses. 

A  wave  energy  spectrum  Is  selected  to  represent  the  random  sea.  Wave 
spectra  are  described  In  Appendix  A.  The  wave  spectrum  represents 
the  distribution  of  the  random  sea's  energy  among  an  Infinite  set  of 
regular  waves  that  when  added  together  create  the  random  character  of 
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the  sea.  By  assuming  that  the  response  is  linear,  the  response  of 
the  offshore  structure  to  a  regular  wave  Is  equal  to  the  RAO  times 
the  regular  wave  amplitude.  By  assuming  that  the  response  to  one 
wave  does  not  affect  the  response  to  another  wave,  the  response  of 
the  offshore  structure  to  a  random  sea  Is  the  sum  of  Its  responses  to 
each  of  the  constituent  regular  waves  In  the  random  sea.  The 
response  Is  therefore  a  collection  of  responses  each  with  a  different 
amplitude,  frequency,  and  phase. 

The  energy  of  each  constituent  wave  1s  proportional  to  the  wave 
amplitude  squared.  The  energy  of  the  response  to  a  constituent  wave 
of  the  random  sea  is  proportional  to  the  response  squared,  or  Is 
proportional  to  the  RAO  squared  times  the  wave  amplitude  squared. 
The  response  energy  may  also  be  represented  by  a  spectrum  from  which 
characteristics  of  the  response  may  be  derived.  From  the  response 
spectrum  characteristics  and  the  wave  spectrum  characteristics,  a 
"transfer  function"  can  be  obtained  which  relates  the  response  and 
wave  characteristics. 

B.2.1  Spectrum  Analysis  Procedure 

The  spectral  analysis  procedure  involves  four  steps:  1)  obtaining 
the  response  amplitude  operators,  2)  multiplying  the  wave  spectrum 
ordinates  by  the  RAOs  squared  to  get  the  response  spectrum,  3) 
calculating  the  response  spectrum  characteristics,  and  4)  using  the 
response  spectrum  characteristics  to  compute  the  random  sea  response 
transfer  function. 

The  RAOs  are  usually  calculated  for  a  discrete  set  of  wave 
frequencies,  and  the  discrete  RAOs  are  then  fit  with  a  curve  to 
produce  a  continuous  function.  The  singular  term  "RAO"  is  used  both 
to  signify  a  single  response  amplitude  to  wave  amplitude  ratio  and  to 
signify  the  continuous  function  through  all  of  the  RAOs.  Any 
response  that  is  linearly  related  (proportional)  to  wave  amplitude 
may  be  reduced  to  an  RAO  function.  Typical  responses  are  motions, 
accelerations,  bending  moments,  shears,  stresses,  etc. 
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Multiplication  of  the  wave  spectrum  ordinates  by  the  RAO  squared  is 
simple.  The  two  underlying  assumptions  are  that  the  response  varies 
linearly  with  wave  amplitude  and  the  assumption  that  the  response  to 
a  wave  of  one  frequency  is  independent  of  the  response  to  waves  of 
other  frequencies. 

Response  spectrum  characteristics  are  taken  from  the  shape  of  the 
spectrum  or  are  calculated  from  the  area  under  the  response  spectrum 
and  the  area  moments  of  the  response  spectrum.  Typical 
characteristics  are  significant  response  amplitude,  maximum  response 
amplitude,  mean  period  of  the  response,  and  peak  period  of  the 
response  spectrum. 

The  random  sea  transfer  function  is  the  ratio  of  a  response  spectrum 
characteristic  to  a  wave  spectrum  characteristic.  A  random  sea 
transfer  function  is  usually  presented  as  a  function  of  the  random 
sea  characteristic  period.  A  typical  transfer  function  might  be  the 
ratio  of  maximum  bending  moment  amplitude  per  unit  significant  wave 
height.  The  transfer  function  is  useful  for  estimating  the  response 
to  another  wave  spectrum  with  similar  form  but  different  amplitude. 

B.2.2  Transfer  Function 


A  transfer  function  converts  input  to  output  for  linear  systems.  A 
transfer  function  is  graphically  represented  in  Figure  B-1.  A 
transfer  function  can  relate  motion  response  to  the  height  of 
incident  waves  directly,  or  a  transfer  function  can  relate  motion 
response  to  wave  force,  or  a  transfer  function  can  relate  member 
stresses  to  wave  or  wind  force. 

For  typical  applications  to  the  design  of  offshore  structures,  the 
input  energy  forms  are  waves,  current  and  wind.  The  desired  output 
forms  are  static  displacements,  dynamic  displacements,  and  member 
stresses. 
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B.2.2.1  Equation  of  Motions 


By  assuming  that  the  motions  are  small  enough  that  the  Inertial, 
damping  and  spring  forces  can  be  summed  linearly,  the  equation  of 
motion  can  be  formulated. 

M*X  +  0**  +  K*X  =  F(x,t) 

where  M  Is  the  mass  matrix  which  Includes  the  structure  mass 
properties  plus  the  hydrodynamic  added  mass  effects, 

0  Is  the  linearized  damping  matrix  which  Includes  the 
viscous  damping,  the  wave  damping,  and  the  structural 
damping  effects, 

K  Is  the  stiffness  matrix  which  Includes  the  waterplane 
spring  properties,  the  restoring  properties  of 
moorings  or  tendons,  and  the  stiffness  properties  of 
the  structure  and  any  foundation, 

X  is  the  system  dlsplacanent  vector, 

t  is  the  system  velocity  vector  »  (dx/dt), 

X  Is  the  system  acceleration  vector,  *  (d^x/dt^),  and 

F  is  the  force  vector  which  may  be  calculated  from 

empirical  methods  such  as  Morrison's  equation  or  from 
diffraction  theory  methods. 

The  equations  of  motion  can  be  solved  with  frequency  domain  or  time 
domain  t  chniques.  The  frequency  domain  solution  Involves  the 
method,  c"  narmonic  analysis  or  the  methods  of  Laplace  and  Fourier 
trans* j  .  The  time  domain  solution  Involves  the  numerical  solution 
by  a  time  step  simulation  of  the  motion. 
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B.2.2.2  Response  Amplitude  Operator 

The  solution  of  the  equations  of  motion  result  In  a  transfer 
function.  The  motion  transfer  function  has  an  In-phase  conponent  and 
an  out-of-phase  component.  The  transfer  function  Is  usually 
represented  In  complex  form, 

X(u)  =  A*{XI(a.)  +  1*X0(o»)l 

or  In  angular  form, 

X(  «  )  =  A*(XI*cos(ut)  +  X0*s1n(a)t)l 

where 

X  Is  the  total  response, 

A  Is  the  wave  height, 

XI  is  the  In-phase  component  of  the  response  for  unit 
wave  height,  and 

XO  is  the  out-of -phase  component  of  the  response  for  unit 
wave  height. 

From  this  equation,  the  response  amplitude  operator  (amplitude  per 
unit  wave  height).  Is  found  to  be, 

RAO  *  SORT  (XI^  +  XO^), 

and  the  phase  of  the  harmonic  response  relative  to  the  wave  Is, 

«  -  ATAN  (XO/XI). 

The  response  can  be  written  In  terms  of  the  RAO  and  phase  as. 
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X(u)  »  A*RA0(a))*COS(wt  +  ♦(«)). 


When  a  spectral  analysis  Is  applied  to  the  transfer  function  the 
wave  amplitudes.  A,  become  a  function  of  wave  frequency,  u,  and 
the  X(u)  Is  replaced  by  the  differential  slice  of  the  response 
power  density  spectrum. 

Sp(«)*du  -  [A(w)*RA0(ai)l^ 


or 


SR(tt)*da»  -  a2(u)*RA02(«) 


or 


SR(u)*du  »  S(«.)*du*RA02(oi) 

Thus,  Sf(u»)  «  S(«)*RAo2(oi) 

The  response  spectrum  S(ui)  Is  therefore  just  the  sea  spectrum  times 
the  RAO  squared. 

For  multiple-degree-of -freedom  systems,  there  Is  coupling  between 
some  of  the  motions,  such  as  pitch  and  heave.  For  example,  to  obtain 
the  motion  or  motion  RAO  for  heave  of  a  point  distant  from  the  center 
of  pitch  rotation,  the  pitch  times  rotation  arm  must  be  added  to  the 
structure  heave.  This  addition  must  be  added  with  proper 

consideration  of  the  relative  phase  angles  of  the  pitch  and  heave 
motions,  and  therefore,  such  addition  must  be  performed  at  the 
regular  wave  analysis  stage.  The  combined  heave  (w/pitch)  RAO  can 
then  be  used  In  a  spectral  analysis  to  obtain  the  heave  spectrum  and 
heave  response  characteristics  at  the  point. 
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B.2.3  Wave  Spectra 


The  wave  spectrum  used  In  the  spectral  analysis  may  be  an  Idealized 
mathematical  spectrum  or  a  set  of  data  points  derived  from  the 
measurement  of  real  waves.  When  a  set  of  data  points  are  used,  a 
linear  or  higher  order  curve  fit  Is  employed  to  create  a  continuous 
function.  Custom  wave  spectra  for  specific  regions  are  often 
provided  as  one  of  the  conventional  Idealized  spectra  with  parameter 
values  selected  to  match  a  set  of  measured  wave  data.  For  areas 
where  there  Is  little  wave  data,  wave  height  characteristics  are 
estimated  from  wind  speed  records  from  the  general  area. 

B.2.3.1  Wave  Slope  Spectra 

For  certain  responses,  particularly  the  angular  motions  of  pitch  and 
roll,  the  RAO  Is  often  presented  as  response  angle  per  unit  wave 
slope  angle.  For  these  cases  the  wave  spectrum  In  amplitude  squared 
must  be  converted  to  a  wave  slope  spectrum.  The  maxlnum  slope  of  any 
constituent  wave  of  the  spectrum  Is  assumed  to  be  small  enough  that 
the  wave  slope  angle  In  radians  Is  approximately  equal  to  the  tangent 
of  the  wave  slope.  The  water  depth  Is  assumed  to  be  deep  enough  (at 
least  one-half  the  longest  wave  length)  that  the  wave  length  Is 
approximately  equal  to: 

(g/2it)*T2  or  2irg/w^. 

By  using  the  Fourier  series  representation  of  the  wave  spectrum, 
selecting  one  constituent  wave,  and  expressing  the  wave  equation  In 
spatial  terms  Instead  of  temporal  terms,  the  wave  slope  Is  derived  as 
follows. 


p 

n«  a*cos(2irx/L)  »  a*cos(xM  /g) 
dn/dx  *  -(aM^/g)*s1n(x<u^/g) 
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Idn/dxl 


IMX  ■ 

Squaring  the  equation  to  get  the  slope  squared, 

[dn/dxl^  «  a^*(MVg^) 

Therefore,  the  wave  spectrum  equation  must  be  multiplied  by  (w^/g^) 
to  obtain  the  slope  spectrum.  The  wave  slope  angle  spectrum  Is  the 
wave  slope  spectrum  converted  to  degrees  squared,  I.e.,  multiplied  by 
(180/»)2. 

B.2.3.2  Wave  Spectra  for  Moving  Vessels 

For  self-propelled  vessels  or  structures  under  tow,  the  forward  speed 
of  the  vessel  or  structure  will  have  an  effect  upon  the  apparent 
frequency  of  the  waves.  The  apparent  frequency  of  the  waves  Is 
usually  referred  to  as  the  encounter  frequency.  For  a  vessel  heading 
Into  the  waves  the  encounter  frequency  Is  higher  than  the  wave 
frequency  seen  by  a  stationary  structure.  For  a  vessel  moving  in  the 
same  direction  as  the  waves,  the  encounter  frequency  Is  less  than  the 
wave  frequency  seen  by  a  fixed  structure,  and  if  the  vessel's  speed 
is  great  enough  It  may  be  overrunning  some  of  the  shorter  waves  which 
win  give  the  appearance  that  these  shorter  waves  are  coming  from 
ahead  Instead  of  from  behind. 

The  encounter  frequency  for  a  regular  wave  Is  given  by  the  following 
relationship. 

2 

(i»g  *  u  +  Vtt*'/g 

where  u  1s  the  wave  frequency  In  radians  per  second  as  seen 
from  a  stationary  observer, 

V  is  the  velocity  component  parallel  to  and  opposite  In 
direction  to  the  wave  direction,  and 
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g  1s  the  acceleration  of  gravity  In  units  conpatible 
with  the  velocity  units. 

The  energy  of,  or  area  under  the  curve  of  the  sea  spectrum  must 
remain  constant. 

/Sg(ug)*dMg  a  /S(w)*dw 

Taking  the  derivative  of  the  encounter  frequency  equation  gives  the 
following. 


d«g  =  [1  +  2Vu/g]*du 

Substituting  the  derivative  Into  the  area  Integral  gives  the 
following. 

/Se(Me)*t^  +  2V(ii/g|*dM  •  /S(ci»)*dcii 

Therefore,  equating  the  Integrands  gives  the  relationship  between  the 
encounter  spectrum  and  the  stationary  sea  spectrum. 

Se(“e)  “  S(m)/(1  +  2Vu/g] 

This  equation  Is  required  to  transform  a  stationary  sea  spectrum  to 
an  encounter  spectrum  for  the  purpose  of  Intergrating  the  responses. 

"•o  • 

However,  If  only  the  response  statistics  are  desired,  and  not  the 
actual  response  spectrum,  then  the  same  substitutions  as  above  can  be 
made. 


Sg  «  S/(l  +  2V«/gl 
dhig  a  (1  +  2V«  /gl*du 
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;re^*Se*do.e= 

Therefore,  the  encounter  frequency  need  only  be  used  to  select  the 
response  amplitude  operator  and  the  Integration  Is  still  over  the 
stationary  ^»*equency,  u. 

I.e.,  rg  *  r(Mg)  ■  r(w  +  Vw^/g) 

B.2.3.3  Short-Crested  Seas 


The  usual  mathematical  representation  of  a  sea  spectrum  Is  one- 
dimensional  with  the  random  waves  traveling  In  a  single  direction 
with  the  crests  and  troughs  of  the  waves  extending  to  Infinity  on 
either  side  of  the  direction  of  wave  travel.  A  one-dimensional 
Irregular  sea  Is  also  referred  to  as  a  long-crested  Irregular  sea. 
In  the  real  ocean  the  waves  tend  to  be  short-crested  due  to  the 
Interaction  of  waves  from  different  directions. 

A  two-dimensional  spectrum  (short-crested  sea)  Is  created  from  a 
standard  one-dimensional  mathematical  spectrum  by  multiplying  the 
spectrum  by  a  "spreading  function."  The  most  commonly  used  spreading 
function  1s  the  "cosine-squared"  function. 

f(<i)  *  (2/w)*cos^4» 

where  is  the  angle  away  from  the  general  wave  heading, 
(-ir/2<*<K/2) 

The  cosine-squared  spreading  function  spreads  the  sea  spectrum  over 
an  angle  +/-  90  degrees  from  the  general  wave  heading. 

To  Incorporate  multi -directional  or  short-crested  Irregular  seas  Into 
a  spectral  analysis,  the  RAOs  for  a  range  of  wave  headings  must  be 
obtained.  A  spectral  analysis  Is  performed  for  each  heading  using 
the  one-d1mens1ona1  sea  spectrum.  The  results  of  the  one-dimensional 
analyses  are  then  multiplied  by  Integration  factors  and  summed. 
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The  following  Is  a  sample  of  a  set  of  heading  angles  and  the 
Integration  factors  for  a  cosine  squared  spreading  function. 


Factor 

0" 

0.2200 

±20* 

0.1945 

±40* 

0.1300 

±60* 

0.0567 

±80* 

0.0088 

B.2.4  Force  Spectrum 

For  simple  single-degree-of -freedom  systems,  a  force  spectrum  can  be 
generated  directly  from  the  calculated  or  measured  regular  wave 
forces. 

The  force  on  the  structure  Is  calculated  by  empirical  or  theoretical 
methods,  or  Is  derived  by  analyzing  measured  strain  records  from 
tests  on  the  structure  or  on  a  model  of  the  structure.  This  force  Is 
the  right  hand  side  of  the  equation  of  motion  as  described  In  Section 
B.2.2.1. 


The  force  itself  has  an  In-phase  and  an  out-of -phase  component 
relative  to  the  regular  wave  which  generates  the  force.  The  force 
can  be  written  In  complex  form, 

F(w)  -  A*lFI(<a)  +  1*F0(«)1 

or  In  force  RAO  and  phase  form, 

F(iii)  »  A*RAO^(ui)*cos(wt  +♦(«)) 

where  RAO^  »  SQRT  (FI^  +  FO^),  and 

«  -  ATAN  (FO/FI). 
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The  force  spectrum  can  be  .created  by  multiplying  a  selected  wave 
spectrum  times  the  force  RAO  squared. 

S^(u))  -  S(m)*RA0^^(«) 

B.2.5  White  Noise  Spectrum 

Most  sea  spectra  have  a  well  defined  peak  of  energy  and  the  energy 
trails  off  to  near  zero  away  from  the  peak.  Other  environmental 
Inputs  that  are  described  by  spectra,  such  as  wind  force,  may  not 
have  a  definite  peak  and  may  even  appear  constant  over  a  wide  range 
(broad  band)  of  frequencies. 

Often  the  response  RAO  Is  narrow  banded,  that  Is,  the  structure  tends 
to  respond  at  a  narrow  range  of  frequencies,  centered  about  a 
resonant  frequency.  When  the  combination  of  a  broad  banded 
excitation  spectrum  and  a  narrow  banded  RAO  exist,  the  spectral 
analysis  can  be  greatly  simplified. 

A  broad  banded  spectrum  can  be  approximated  by  a  "white  noise 
spectrum"  which  has  constant  energy  over  the  whole  frequency  range  of 
the  spectrum. 

For  a  single  degree  of  freedom  system,  the  response  can  be  defined  In 
terms  of  a  "dynamic  amplification  function"  times  an  expected  static 
displacement.  The  dynamic  amplification  function  Is  as  follows, 

IH(ai)|  -  l/((W»n)')* 

where 


u  Is  radial  frequency, 

u  1s  the  undamped  "natural  frequency", 

*  (k/m)**. 
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Z  Is  the  damping  ratio*  the  ratio  of  the  actual  damping  to 
the  critical  damping.  ^  >  c/(4lcm)^, 

k  Is  the  spring  constant, 

m  Is  the  mass  that  Is  In  motion,  and 

c  Is  the  actual  damping. 

The  expected  static  displacement  Is  simply  force  divided  by  the 
spring  constant,  or  the  expected  static  displacement  spectrum  Is  as 
follows, 

S,(.)  -  S,(.)/k’ 

From  these  equations,  the  response  spectrum  Is  found  to  be, 

R(.)  -  (l/k’)*IH(«)l**S,(»), 
and  the  mean  squared  response  Is, 


y’  (t)  =,r(l/k)  **  IH(«)  I 

The  (l/k)‘  Is  constant,  and  by  approximating  the  force  spectra,  by  a 
white  noise  spectrum  with  magnitude  Sf(u^),the  mean  squared  response 
Is  simplified  to. 


y'(t)  -  (s,(«)„/k’)*^r  IH<.)l’*d.. 

For  lightly  damped  systems,  (^«1),  the  Integral  may  be  evaluated  to 
yield. 


y'(t)  •  l.*.„*Sf(.„)l/(2*t*k') 
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B.3 


EXTREME  RESPONSE 


The  extreme  response  of  an  offshore  structure  may  be  determined  In 
two  ways.  An  extreme  environmental  event  may  be  selected,  and  the 
responses  to  the  extreme  event  then  calculated.  A  set  of 
environmental  spectra  can  be  selected;  the  response  spectra  to  each 
environmental  spectra  calculated;  and  the  extreme  responses  derived 
by  statistical  analysis  of  the  response  spectra.  The  first  method  Is 
often  called  a  "deterministic*  method,  and  the  second  method  Is 
referred  to  as  a  "probabilistic"  method.  In  actual  design  practice 
the  two  methods  are  often  Intermixed  or  cnnblned  In  order  to  confirm 
that  the  extreme  response  has  been  found. 

B.3.1  Maximum  Wave  Height  Method 

In  deterministic  design,  a  set  of  extreme  conditions  Is  supplied  by 
oceanographers  or  meterologists.  The  extreme  conditions  are  of 
course  derived  from  statistical  analyses  of  wave  and  weather  records, 
but  the  design  engineer  Is  usually  not  Involved  In  that  stage  of  the 
calculations. 

The  given  extreme  conditions  are  applied  to  the  offshore  structure  to 
determine  the  various  responses.  Unfortunately,  the  given  extreme 
conditions  may  not  always  produce  the  extreme  responses.  For 
example,  the  prying  and  racking  loads  governing  the  design  of  many 
structural  members  of  semisubmersibles  are  typically  maximized  In 
waves  with  lower  heights  and  shorter  lengths  then  the  maximum  height 
wave.  Tendon  loads  on  tension  leg  platforms  (TLPs)  are  also  often 
maximized  In  waves  that  are  lower  and  shorter  than  the  maximum  wave. 

Since  the  oceanographer  or  meterologist  who  produced  the  set  of 
extreme  conditions  does  not  have  Information  about  the 
characteristics  of  the  offshore  structure,  he/she  Is  unable  to  select 
an  extreme  or  near  extreme  condition  that  will  produce  the  greatest 
response.  Conversely,  the  design  engineer  usually  has  little  or  no 
Information  about  the  wave  and  weather  data  that  was  used  to  derive 


B-15 


the  set  of  extreme  conditions,  and  thus,  he/she  Is  unable  to  create 
alternate  conditions  to  check  for  greater  response. 

The  design  engineer  may  request  a  range  of  extreme  conditions,  such 
as:  the  maximum  height  wave  with  a  period  of  9  sec,  the  maximum 
height  wave  with  a  period  of  10  sec,  etc.  The  Increased  number  of 
conditions  Increases  the  number  of  analyses  required,  but  allows  the 
design  engineer  to  confirm  which  conditions  produce  the  extreme 
responses. 

The  maximum  wave  height  method  Is  best  used  when  the  response  Is 
highly  nonlinear  and  the  spectral  analysis  method  Is  therefore  not 
appropriate. 

B.3.2  Wave  Spectrum  Method 

A  full  probabilistic  analysis  Involves  calculating  responses  to  the 
entire  suite  of  possible  environmental  conditions.  Statistical 
analysis  of  these  responses  Is  then  performed  In  order  to  predict  a 
suitable  extreme  for  each  response.  This  requires  far  fewer 
assumptions  on  the  part  of  those  who  supply  environmental  criteria, 
but  a  much  more  extensive  set  of  environmental  data. 

With  the  wave  spectrum  method,  a  set  of  wave  spectra  are  provided  by 
oceanographers  or  meterologists.  The  RAOs  for  the  response  of 
Interest  are  squared  and  multiplied  by  the  wave  spectrum.  A  wave 
spectrum  Is  assumed  to  represent  a  Gaussian  random  distribution. 
Since  the  response  spectrum  Is  created  by  a  linear  multiplication, 
the  response  spectrum  also  represents  a  Gaussian  random 
distribution.  The  significant  response,  maximum  response,  etc.  can 
be  calculated  using  the  equations  for  calculating  the  significant, 
maximum,  etc.  wave  heights. 

The  equations  for  maximum  wave  height  are  summarized  here  In  terms  of 
response: 
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significant  response,  (DA): 


Rj  =  4.00*(nio)** 

Maximum  response  In  1000  cycles,  (DA): 

*^1/1000  “  7.43*(mQ)^  ■  1.86*R5 

Maximum  response  Is  t  hours,  (DA): 

W  ’  2*(2*mo*ln(3600n/T2)]*» 

where  mg  Is  the  area  under  the  response  spectrum, 

T2  Is  the  zero-up-crossIng  period  of  the  response  as  found 

from  the  equation, 

Tj  »  2ir*(mg/m2)^,  and 

m2  Is  the  second  radial  frequency  moment  of  the  response  spectrum. 
B.4  OPERATIONAL  RESPONSE 


In  order  to  determine  the  normal  day-to-day  motions  and  stresses  to 
assess  motion  related  downtime  and  fatigue  damage,  the  distribution 
of  wave  heights  versus  wave  periods  are  considered.  A  wave  scatter 
diagram  condenses  and  summarizes  wave  height  and  wave  period 
statistics.  It  Is  a  two-parameter  probability  density  function. 
Typically  a  wave  scatter  dlagrwa  Is  presented  as  a  grid  of  boxes, 
with  one  axis  of  the  grid  being  average  zero-up-cross  1ng  periods  and 
the  other  axis  being  significant  wave  heights.  Within  the  boxes  of 
the  wave  scatter  diagram  are  numbers  which  represent  the  percentage 
of  the  sea  records  having  the  corresponding  characteristics  of  Hs  and 
Tz  see  Figure  A-2. 
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A  response  scoter  diagram  could  be  generated  by  taking  the  wave 
spectrum. for  each  sample  used  to  create  the  wave  scatter  diagram  and 
performing  a  spectral  analysis  for  the  response.  The  computed 
characteristics  are  then  used  to  assign  the  percentage  of  occurrence 
to  the  appropriate  box  In  the  response  scatter  diagram.  This  entails 
considerable  work  and  can  be  simplified  by  reducing  the  number  of  sea 
spectra  considered,  as  described  below. 

B.4.1  Special  Family  Method 

All  of  the  original  sea  spectra  used  to  define  the  wave  scatter 
diagram  must  be  available.  In  order  to  select  a  special  family  of  sea 
spectra  to  represent  the  whole  population. 

The  sea  spectra  are  first  grouped  by  wave  height  bands,  such  as  0  to 
2  ft  significant  wave  height,  2  ft  to  4  ft  H^,  etc.  The  average 
properties  of  the  spectra  within  a  group  are  computed.  Within  each 
group,  which  may  contain  thousands  of  sample  sea  spectra,  a  small  set 
of  sea  spectra  are  selected  to  represent  all  of  the  spectra  In  the 
group.  The  small  set  will  typically  contain  4  to  10  spectra. 

The  spectra  of  a  representative  set  are  selected  by  a  Monte  Carlo 
(Shotgun)  process  which  randomly  picks,  say  8,  spectra  from  the 
group.  The  mean  spectrum  and  the  standard  deviation  of  the  spectral 
ordinates  about  the  mean  spectrum  are  computed  for  the  8  spectra.  A 
weighted  sum  of  differences  In  properties  between  the  8  spectra  and 
the  total  population  of  the  group  represent  the  "goodness  of  fit"  of 
that  set  of  8  spectra. 

A  second  representative  set  of  8  spectra  is  then  selected  from  the 
group,  and  the  "goodness  of  fit"  of  the  second  set  Is  computed.  The 
better  set  (first  or  second)  Is  retained  and  compared  to  a  third 
sample  of  8,  etc.  The  process  Is  repeated  many  times,  say  1000, 
within  each  wave  group. 


8-18 


From  this  process,  the  original  number  of  sea  spectra,  which  may  have 
been  thousands.  Is  reduced  to  the  number  of  wave  height  bands  times 
the  number  of  spectra  In  each  representative  set. 

B.4.2  Wave  Spectrum  Method 

A  reduced  set  of  sea  spectra  can  be  generated  to  represent  the 
variation  of  and  T2  as  given  In  a  wave  scatter  diagram. 

If  the  original  sea  spectra  are  not  available,  a  set  of  sea  spectra 
can  be  created  directly  from  the  wave  scatter  diagram.  In  this  case 
the  shape  of  the  spectrum  must  be  assumed.  For  various  areas  of  the 
world's  oceans,  preferred  mathematical  spectrum  equations  exist. .  For 
the  North  Sea,  the  mean  JONSWAP  spectrum  Is  preferred.  For  open 
ocean,  the  Bretschnelder  (ISSC)  spectrum  Is  preferred.  For  the  Gulf 
of  Mexico,  the  Scott  spectrum  has  been  recommended. 

Using  the  and  for  each  populated  box  In  the  wave  scatter 
diagram,  and  the  selected  sea  spectrum  equation,  a  set  of  wave 
spectra  are  defined.  With  this  method  the  number  of  sea  spectra  Is 
reduced  to  the  number  of  populated  boxes  In  the  wave  scatter  diagram, 
but  no  more  than  the  number  of  wave  height  bands  times  the  number  of 
wave  period  band. 
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STRESS  CONCENTRATION  FACTORS 


C.l  OVERVIEW 

C.1.1  Objectives  and  Scope 

A  comprehensive  document  on  stress  Concentration  factors  (SCF)  would 
Include  assessment  of  test  results,  detailed  review  of  empirical 
equations,  evaluation  of  finite  element  studies,  and  presentation  of 
parametric  studies  showing  the  sensitivities  of  parameters  affecting 
SCFs. 

The  objective  of  this  appendix  Is  limited.  Following  a  brief 
discussion  of  empirical  equations,  parametric  study  results  are 
presented  to  assist  the  engineer  1n  avoiding  undesirable  Joint 
details.  The  sensitivity  and  Interaction  of  variables  shown  In 
tables  and  figures  also  allow  quick  assessment  of  steps  necessary  to 
Improve  other  geometries. 

Empirical  formulations  are  applicable  to  a  limited  range  of  simple 
Joint  geometries.  A  complex  Joint  often  requires  carrying  out  of  a 
finite  element  analyses  (FEA)  to  determine  the  SCFs.  The  results  of 
a  FEA  Is  also  presented  to  Illustrate  the  applicable  SCFs  for  a  given 
geometry. 

Cl. 2  Current  Technology 

The  SCF  values  can  be  computed  through  the  use  of  a  number  of 
alternative  equations.  These  equations  have  been  mostly  based  on 
analytical  (finite  element)  and  small-scale  experimental  (acrylic 
model  test)  work.  The  tests  carried  out  on  Joints  that  reflect  those 
In-service  (I.e.  both  In  size  and  fabrication  methods)  are  few  and 
limited  to  several  simple  Joint  configurations.  Thus,  the  equations 
available  should  be  reviewed  carefully  to  ascertain  their  range  of 
validity  and  overall  reliability  prior  to  their  use  In  design. 
Considering  the  simple  Joint  configurations  of  T,  Y,  OT,  K  and  X,  the 
equations  available  for  use  In  design  are: 


C-1 


0  Kuang  (Reference  C.l) 

0  Wordsworth  (References  C.2,  C.3) 

0  Gibstein  (References  C.4,  C.5) 

0  Efthymlou  (Reference  C.6) 

0  Marshall  (Reference  C.7) 

0  UEG  (Reference  C.8) 

There  are  significant  differences  In  the  validity  ranges  of  these 
equations.  The  SCFs  computed  based  on  different  equations  also  often 
vary  considerably.  The  Kuang  equations  are  applicable  to  T,  Y,  and  K 
joints  for  various  load  types.  Wordsworth  and  Woodsworth/Smedley 
equations  are  applicable  to  all  simple  Joints.  Gibstein  equations 
are  applicable  to  T  joints  while  the  Efthymlou  equations  cover  T/Y 
joints  and  slmple/overlapplng  K/YT  joints.  The  equations  proposed  by 
Marshall  are  applicable  to  simple  joints,  based  on  those  equations  by 
Kellogg  (Reference  C.9),  and  were  Incorporated  Into  API  RP  2A. 

Substantial  work  has  been  carried  out  to  validate  the  applicability 
of  various  SCF  equations.  Although  some  of  the  work  carried  out  by 
major  oil  companies  are  unpublished,  such  work  still  Influence  on¬ 
going  analytical  and  experimental  research.  Delft  von  O.R.V.  et  a1. 
(Reference  C»10)  indicate  that  the  UEG  equations  offer  a  good 
combination  of  accuracy  and  conservatism  while  the  Efthymlou  (I.e., 
She11-SIPM)  equations  show  a  good  comparison  with  experimental 
data. 

Ma  and  Tebbet  (Reference  C.ll)  report  that  there  Is  no  consensus  on 
whether  a  design  SCF  should  represent  a  mean,  lower  bound  or  some 
other  level  of  confidence.  Tebbett  and  Lalani's  (Reference  C.12) 
work  on  reliability  aspects  of  SCF  equations  Indicate  that  SCF 
equations  underpredicting  the  SCF  values  In  less  than  16X  of  the 
cases  can  be  considered  reliable.  Thus,  when  presenting  the  findings 
of  45  elastic  tests  carried  out  on  15  tubular  joints  representing 
typical  construction,  Ma  and  Tebbet  report  that  Wordsworth,  UEG  and 
Efthymlou  equations  meet  this  criteria  and  offer  the  best 
reliability. 
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Ma  and  Tebbett  also  state  that  while  both  UEG  and  Wordsworth 
equations  overpredict  X  joint  SCFs,  none  of  the  equations  overpredict 
the  K  Joint  SCFs.  The  comparative  data  Indicate  that  the  SCFs 
computed  using  Kuang  and  Gibstein  equations  for  T/Y  Joints  subjected 
to  axial  loading  under  predict  the  measured  data  In  more  than  16X  of 
the  cases.  (See  Figure  C.1-1). 

Tolloczko  and  Lalani  (Reference  C.13)  have  reviewed  all  available  new 
test  data  and  conclude  that  reliability  trends  described  earlier  for 
simple  Joints  remain  valid  and  also  state  that  Efthymiou  equations 
accurately  predict  the  SCFs  for  overlapping  Joints. 
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C.2 


STRESS  CONCENTRATION  FACTOR  EQUATIONS 


C.2.1  Kuanq  with  Marshall  Reduction 

The  Kuang  stress  Concentration  factor  equations  for  s1^)1e 
unstiffened  Joints  are  shown  on  the  following  page.  The  brace  stress 
Concentration  factor  equations  Include  Marshall  reduction  factor, 
o'*.  The  validity  ranges  for  the  Kuang  stress  Concentration  factor 
equations  are: 


Tens  Validity  Range 


d/0 

0.13  -  1.0 

T/D 

0.015  -  0.06 

t/T 

0.20  -  0.80 

g/0 

0.04  -  1.0 

D/L 

0.05  -  0.3 

e 

25  >  90 

where,  0  •  chord  dlaaeter 
T  •  chord  thickness 
d  -  brace  dlaaeter 
t  «  brace  thickness 
g  >  gap  between  adjacent  braces 
L  ■  chord  length 

e  >  angle  between  brace  and  chord 


C-4 


C.2.2  Smed 1 ey-Wordswort h 

The  Smedley-WordSMorth  stress  Concentration  factor  equations  for 
simple  unstiffened  Joints  are  shown  on  the  following  pages.  The 
notes  for  the  equations  shown  on  the  following  pages  Include  the 
Shell  d/0  limitation  of  0.95.  This  Interpretation  Is  open  to  a 
proJect-by-proJect  review. 

The  validity  ranges  for  the  Smedley-Uordsworth  equations  are: 


Term  Validity  Range 


d/0 

0.13  -  1.0 

0/2T 

12.0  -  32.0 

t/T 

0.25  -  1.0 

9/0 

0.05  -  1.0 

2L/0 

8.0  -  40 

30  >  90 

where.  0  •  chord  diameter 
T  »  chord  thickness 
d  ■  brace  diameter 
t  «  brace  thickness 
g  «  gap  between  adjacent  braces 
L  ■  chord  length 
«  angle  between  brace  and  chord 
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PARAMETRIC  STUDY  RESULTS 


Figures 

The  Kuang  and  Snedley-Mordsworth  chord  stress  Concentration  factors 
for  T  Joints  are  shown  In  Section  C. 3. 1(a)  and  C. 3. 1(b), 
respectively.  The  Kuang  and  Smedley-Wordsworth  chord  stress 
Concentration  factors  for  K  Joints  are  shown  In  Section  C. 3. 1(c)  and 
C. 3. 1(d),  respectively.  The  Saedley-Mordsworth  chord  stress 
Concentration  factors  for  X  Joints  are  shown  In  Section  C. 3. 1(e). 
Since  the  chord  side  of  the  weld  stress  Concentration  factor  1s 
generally  higher  than  the  brace  side  of  the  weld  stress  Concentration 
factor,  only  the  chord  side  of  the  weld  stress  Concentration  factors 
are  shown. 


C.3.1(a)  Kuang  Chord  SCF’s  for  T-Joints 


The  Kuang  chord  SCF's  for  T-joints  are  shown  on  the  following 
pages.  The  following  parameters  are  assumed  for  the  Kuang  figures: 

1)  Y  »  D/2T  »  12.0 

2)  Q  •  «  30.0  degrees 

3)  -  *  0/L  »  0.0571 
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Out  aV  fUnt  icaaiay  $Cf 


Chord  Side 


Bract  Sidt 


K^Joints 

SCF„  . 

SCFgy  -  1.4007‘°-®/!0-0®T°'®^slnO-®9 

T-Joints 

SCF„  -  1.177T-°-*»e-1-2<^Tl-333a-0J)67,i„1.»4g 
SCFgy  -  0.4637'°  ®|J’°'°*T°-®®slnO-®^d 
SCFgj  -  0^77*’  0^^T0-®®®^0-^®^$in 

forOJ<d<0^ 

SCFgj  -  0^97  ’  °^^r0*®®®/3'0‘®^®iln 
for  0jS</3  <0.75 


SCFb,  -  0.8257-0-^®’^‘0-^’r0-5»i,0-«8,^-^«^ 
SCF^  -  2.827/J'°-*T0**iln0-®fl 

SCF,j^  -  2.7847‘®-®®s'’-®^  ra^*'2  jj„  1  J4  j 

SCF^  -  1.1097‘°‘^|J'°-“rO*®ilnO-2'5 

SCF^^  -  O.8437^^rO•5«0O•«>^in2J^33^ 

for0J</3<0£5 

SCFb,  -  0.441  7'O-®®2rO'5«0-O-281, in 2J)33^ 
for  0£5<|3  <0.75 


TK-Joints 

SCFjj^  -  1.287^^^®‘‘2r^-®®®iind 


SCF^^  -  6.6B7'°’'^  w0-^2e^„0J9 

forO*<B<46’ 

SCFj,^  -  1238  7'°-'^'°'2®t0  ®  wO*^26jjn2J8  5 


®^^bx(c»mr»l) 


for45’<e<90* 

4.4918  7  •0  '‘23p-0J96  ^  0.672  ^^0.159  2.257  q 


Whtr* 

7  ■  chord  thiekntn/chord  diarnttar 

6  *  angit  bttwatn  braoa  and  chord 

r  >  braoa  thicknan/chord  thicknats 

a  ■  chord  diamatar/chord  langdi  batwaan  tupporta 

V  *  saparation  distanca  batwaan  bracas/chord  diamatar 

CP  ■  saparation  distanca  batsraan  bracas  for  TK  Joints/chord  diamatar 

0  •  braoa  diamatar/chord  diamatar 


Kuang  Formulas  Usad  for  Computing  Stress  Concentration  Factors 


Chord  Side 


Brace  Side 


K-Joints 

SCFgjj  -  1.8  (Tsin^v'r  ) 

SCFbx  -  1 .0  +  0.6  Or  (1 .0  +v/|.  SCF^jl  >  1  -8 

SCFgy  »  1 .2  (t  sin  0  ) 

SCFby-  1.0  +  0.6Qr[1.0+>/|.SCFcyl  >1.8 

SCFgj  *  2.7  (t  sin  6  ) 

SCFbj  -  1.0  +  0.6Qr  ll.0+>/|.SCFcj]  >1.8 

Y-Branch  Joints 

SCF|^uang“  2.06  7°  ®°®  (sin0)^  ®°^  r^-333 

SCF^yy^  "  for  7  <25 

■  1 .5  r  sin  0  7°*^  for  7  >  25 

SCFcxmod-  SCFc,  +  rcos0 

SCFjc  -  SCF|^u3ng  < 

SCF-n,-  1.0 +  0.6  Or  I1.0+v/^.SCFjgl  >1.8 

'^SCFgjj^jyj 

>SCFb^ 

SCFy  »  same  as  for  K 

SCFy  ■  same  as  for  K 

SCFj  »  same  as  for  K 

SCFj  *  same  as  for  K 

Unreinforced  Cross  Joints 

SCF„  -  1.333  ISCFtc)  ^  ^  + 'c 

*  y-branch  ~ 

SCFy  -  1.333  ISCFgy) 

SCFbx-  1.0  +  0.6  0^  (1.0+v/|.SCFjj]  >1.8 
SCFby-  1.0 +  0.6  Or  (1.0 +y|.  SCFy]  >1.8 

SCFj  »  1.333  (SCF„) 

SCFbj  «  1 .0  +  0.6  ( 1 .0  +>/| .  SCFj]  >  1.8 

Where 

B  "  angle  between  brace  and  chord 
0  "  can  diameter 
T  ■  can  thickness 
tg  ■  nominal  chord  thickness 
d  *  brace  diameter 
t  •  stub  thickness 
r  «  t/T 
7  -  (D-T)/2T 
0  -  d/D 

-  exp  (-(0.5  T  +  t)/5s/5Tdtl 


Marshall  Formulas  Used  for  Computing  Stress  Concentration  Factors 


Chord  Sidt 

SCP^  -  1.77t^(2.42  -  2^  -  14.4/J)^ 

SCFgy  -  0.757°*®t0*8(1.6^^/^  -  0.7^2)sln(i.5  - 
SCFgj-TTPd.Se  -  1.46^)8ln^^<^5  -  14.4/J)^ 

BnceSidft 

SCFbx-1  +  0.63SCFejj 

SCFby-l  +  0.63SCFgy 
SCFbj-1  +  0.63SCFgj 


When 


p  ■  Brace  Olameter/Chord  Diameter 
7  ■  Chord  Radiuf/Chord  Thicknefi 
r  ■  Brace  Thickness/Chord  Thicknes 
S  -  Acute  Angie  Between  Brace  and  Chord 


Smediey  Formula*  Used  for  Computing  Stress  Concentration  Factors 
for  Unninforced  Cross  Joints 


Defioicion  of  Parameters.  Validicv  Ranges  and  Notes  on  Tables 


Oefinicion  of  Tubular  Joinc  Paraaccers 


~  2L/0  where  D  ■  chord  outside  dienecer 
P  >  d/D  I  ■  chord  wall  thickness 

s'  >  D/2T  L  •  chord  length  (distance  between  points  of 

contraf lexure) 

T  «  c/T  d  ^  brace  outside  diameter 

t  ■  brace  wall  thickness 

^  >  g/0  g  ■  gap  between  adjacent  braces 

Validity  Ranees  for  Parametric  Equations 

a  <  -^  <  40 
0.13  <  p  <  l.O 
12  <  y  <  32 

0.25  <  T  <  1.0 

0.01  <  ?  <  1.0  (Tables  I  and  2)  -0.4  <  ?  <  0.01  (Table  3) 

30*  I  8  <  90* 

0,  >  07 

63  <  90* 

N’otes  on  Tables 


(1)  Tables  1  to  3 

(1)  SCF  (brace)  •  1  4-  0.63  SCF  (chord) 

(2)  SCFs  are  lioiiced  cu  a  minimum  value  of  1.6 

(3)  For  joints  outside  specified  validity  ranges  calcuiats  SCF 
actual  joint  parameters  and  with  parameter*? 

validity  limit.  The  greater  SCF  value  obtained  should  be  used 
in  analysis.  See  also  notes  B(l)  and  C(l)  below. 

(2)  Table  I  only 

(1)  If  0.98  £  P  5  1*0  then  use  P  ■  0.98. 

i'i)  The  K  and  XT  joint  equations  are  based  on  nominal  stress  in 
brace  1 . 

(3)  For  XT  joints  where  the  load  in  brace  3  is  smaller  than  10*  of 
the  maximum  load  in  adjacent  braces  I  and  2  the  joL:iL  t/pe 
should  be  re-categorised  as  X  with  g  the  gap  oecveen  oraces  i 
and  2 . 

(4)  The  equations  indicated  for  K  and  XT  joints  apply  only  to 
balanced  axial  load. 


Tabic  2  only 

(1)  If  P  >  0.95  for  oue-of-planc  bending  chen  use  p  -  0.95. 

(2)  ^e  equations  indicated  for  K  and  KT  joints  apply  only  ta 

bMdinl.**"  **‘^****  ouc-of-plane 

Table  3  only 

?;!r°f:nr.s:rri 


TYPE  AND  LOAOINi}  JiUCRO.  SADDLE  SCF  CHORD.  CROWN  SCF 


Pb(i>mel'ic  Equalicni  lor  SCF'»  cl  Ur.}iill«(ic<i,  Nori-owciloppctl  Tuboior  JainU  •  Aiiol  Loariinq  in  Broct 


Kuong  SCF  Computation 


Kuong  SCF  Computation 


C. 3. 1(b)  Sraedley-Wordsworth  Chord  SCF's  for  T-Joints 

The  Smedley-Wordsworth  chord  SCF's  for  T-joints  are  shown  on  the 
following  pages.  The  following  parameters  are  assumed  for  the 
Sraedley-Wordsworth  figures: 

1)  Y  =  0/2T  -  12.0 

2)  G  »  «  30.0  degrees 

3)  .  =  2L/D  *  35.0 

The  Shell  d/D  limitations  have  not  been  Imposed  for  the  SCF 
calculation. 
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Smedley-Wordsworth  SCF  Computation 


Smediey-Wordsworih  SCF  Computation 


Rlan«  SCR  Crown  Position 


Smediey-Wordsworth  SCF  Computation 
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Smedley-Wordsworth  SCF  Computation 


C. 3. 1(c)  Kuang  Chord  SCF's  for  T-Joints 

The  Kuang  chord  SCF's  for  K- joints  are  shown  on  the  following 


pages.  The  fol 

lowing  parameters 

are  assumed  for  the  Kuang  figures: 

1) 

Y  =  D/2T  * 

12.0 

2) 

0  =  » 

30.0  degrees 

3) 

«  =  D/L  = 

0.0571 
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Axtol  SCF 
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Kuang  SCF  Computab'on 


Kuang  SCF  Computation 


Kuong  SCF  Computation 


C. 3. 1(d)  Smedley-Wordsworth  Chord  SCF's  for  K-Joints 

The  Smedley-Wordsworth  chord  SCF's  for  K- Joints  are  shown  on  the 
following  pages.  The  following  parameters  are  assumed  for  the 
Smedley-Wordsworth  figures: 

1)  Y  -  0/2T  =  12.0 

2)  0  »  ®2  “  degrees 

3)  «  *  2L/D  «  35.0 

The  Shell  d/D  limitations  have  not  been  Imposed  for  the  SCF 
calculation. 
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Rlan«  SCF - Crown  Position 
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Smedley-Wordsworth  SCF  Computation 
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Smediey-Wordsworth  SCF  Computation 
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C. 3. 1(e)  Sraedley-Wordsworth  Chord  SCF's  for  X-Joints 

The  Smedley-Wordsworth  chord  SCF's  for  X- Joints  are  shown  on  the 
following  pages.  The  following  parameters  are  assumed  for  the 
Smedley-Wordsworth  figures: 

1)  Y  =  D/2T  =  12.0 

2)  0  »  =  30.0  degrees 

3)  .  =  0/L  =  35.0 

The  Shell  d/D  limitations  have  not  been  imposed  for  the  SCF 
calculation. 
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C.3.2  Tables 


The  Kuang  and  Smed ley- Wordsworth  chord  stress  Concentration  factors 
for  T  joints  are  shown  in  Section  C. 3. 2(a)  and  C. 3. 2(b)* 
respectively.  Since  the  chord  side  of  the  weld  stress  Concentration 
factor  is  generally  higher  than  the  brace  side  of  the  weld  stress 
Concentration  factor*  only  the  chord  side  of  the  weld  stress 
Concentration  factors  are  shown. 


C-12 


C.3.2(a)  Kuang  Chord  SCF's  for  T-Joints 

The  Kuang  chord  SCF's  for  T-joints  are  shown  on  the  following 
pages.  The  following  parameters  are  assumed  for  the  Kuang  figures: 

1)  -  =  D/L  =  0.0571 
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T-Jriik  Mai  BF 
CM  Sidi  of  Wi 


6mi>  12.0  6Ma>  110  GHm*  20.0  BMi>.2S.O 

j - 1 - 1 - 1 - 1 - 1 

1Ma>iTM*S  MtidA  Ma«M  I  Ma^d/I  Ma4/I  1 

I  t/r  t  0.3  1  0.S  I  0.7  !  0.9  I  0.3  1  0.S  I  0.7  1  0.9  I  0.3  1  0.S  I  0.7  1  0.9  I  0.3  1  0.S  1  0.7  1  0.9  I 

■  ■......*.■■•11  I _ I _ I _ I 


i  0.20  10.032  0.302  0.432  0.272  10.737  0.073  0.318  0.320  10.133  0.819  0.031  0.411  11.144  1.017  0.783  0.192  I 

1 - 1  1111 

30.0  dql  0.40  11.393  1.410  1.090  0.080  11.908  1.010  1.300  0.821  12.407  2.140  1.048  1.037  12.883  2.SB  1.973  1.241  1 

0.324  ridl - 1  1111 

1  0.00  12.730  2.432  1.872  1.178  13.270  1913  1212  1.411  14.134  3.073  1829  1.780  14.130  4.401  3.388  1132  1 

I - 1  I  1  1  1 

I  0.80  14.014  3.309  1747  1.729  14.808  4.274  3.290  1070  10.000  1395  4.131  1012  17.204  0.438  4.972  113  1 


I  0.3  11.137  1.011  0.778  0.490  11.302  1.211  0.9S?  0.380  11.719  !.S28  1.170  0.740  11038  1.830  1.409  0.S  I 

! - !  1111 

43.0  dqi  0.40  11800  1348  1.901  1.234  13.432  3.082  1399  1.478  14.331  3.830  1104  1.803  13.187  4.011  3.330  1233  1 

0.783  rail - !  1111 

1  0.00  14.921  4.373  3.308  1119  11895  1239  4.085  13B  17.430  0.011  1089  1202  11903  7.917  0.094  1833  1 

1 - 1  1111 

1  180  17.221  0.420  4.912  1110  IIOIB  7.088  1919  1721  110.91  9.700  7.407  4.099  11100  11.01  1943  1027  1 


1  120  11.004  1.420  1.097  0.090  11.921  1.707  1.314  1B7  12.423  2.134  l.OSB  1.043  12.902  2.380  1.180  1.230  1 

I  1 - 1  1111 

I  00.0  dagi  0.40  14.041  1312  2.705  1.740  14.851  4.302  1312  2.084  10.100  1428  4.179  2.029  17.312  0.301  1004  1149  1 

!  1.047  raril - 1  1  1  1  1 

!  0.00  10.938  0.108  4.748  2.988  10.308  7.387  1080  1378  110.48  9.320  7.174  4.314  112.35  11.10  LS9Z  1400  1 

!  1 - 1  1111 

I  0.80  110.18  9.051  0.907  4.384  112.19  10.83  8.3M  120  11138  1107  10.32  0.02  11142  10.37  12.00  7.933  1 


!  12  12.040  1.819  1.400  1881  12.431  2.179  1.077  1.035  11092  2.749  LUO  1.331  11703  1212  LS34  1.315  1 

! - 1  1111 

90.0  dagt  140  11130  4.384  132  L22  10.12  1481  4.220  L29  17.790  0.92  1332  1335  19.32  1215  0.385  4.018  i 

1.571  radl - !  1111 

!  100  11832  7.82  0.038  1812  11100  9.42  7.25  4.305  11137  11.89  9.154  1700  110.01  14.2  1191  0.818  1 


Ite|  8f  Dapatatin 


T-Jsiak  b-PlMi  SCF 
Qvri  Sidi  of  Mold 


I  bai>  12.0  GMi>  110  i  6mh>  20.0  !  tei>  2.0  I 

! - 1 - 1 - 1 - 1 - 1 

1hrta«!Tai>i  Miid/I  Ma-ri/0  i  Mo>ri/D  I  Bili^ 

!  t/r  I  0.3  i  0.S  I  0.7 :  0.9 1  0.3  s  o.s  i  0.7 1  0.9 1  0.3 :  o.s  \  0.7 :  0.9 1  0.3 1  o.s  i  0.7 :  0.9 1 

- 1 - 1 - 1 - 1 - 1 — I — I - 1 - 1 - 1 - 1 - 1 - 1 — I — I - 1 — I — I 

I  0.2  I0.3S1  0.SM  0.SS  0.32  10.031  0.011  0.010  O.OOS  :0.749  0.731  0.721  0.717  10.837  0.M0  0.08  0.82  1 

! - 1  1111 

2.0  d0|t  0.40  11.001  0.90  0.908  0.938  11.145  1.122  1.107  1.090  11.301  1.333  1.32  1.308  11.350  1.324  l.SOI  1.4B9  1 

0.32  ridl - 1  1  1  1  1 

1  0.00  11.419  1.391  1.372  1.338  11.02  1.390  1.319  l.SBS  11.929  1.890  1.804  1.840  12.203  2.100  2.131  1110  I 

! - 1  111! 

:  0.2  11.818  1.781  1.737  1.72  1102  1(«  1009  1.9B9  1142  142  1388  1304  1182  1707  172  1703  1 

I  0.2  10.072  0.0a  0.02  0.043  10.72  0.733  0.743  0.738  10.913  0.895  O.BD  0.82  11.044  1.083  1.009  0.999  I 

■  I - 1  I  1  1  I 

1  4S.0  dqi  0.2  11.22  1.195  1.12  1.12  11.32  1.307  1.32  1.335  11.038  1.02  1.003  1.387  11.890  1.837  1.832  1.814  1 

I  0.7B  r21— — -I  lilt 

1  0.2  11.729  1.094  1.072  1.08  11.977  1.937  1.911  1.892  12.32  2.308  2.272  2.22  12.27  2.0S  2.397  2.Sn  I 

1  5 - 1  I  I  1  I 

■  _ 1  0.2  11213  112  2.141  2.12  12.532  2.42  2.42  2.12  13.20  2.949  2.909  2.82  I3.MI  3.371  3.32  3.213  1 

I  0.2  10.72  0.739  0.729  0.722  10.813  0.845  0.831  0.823  11.085  1.004  0.991  0.92  11.172  1.1«  1.133  1.122  1 

! - 1  1111 

00.0  dqt  0.40  11.370  1.342  1.32  1.311  11.300  1.32  1.24  1.499  11.801  1.8S  1.799  1.72  12.12  2.08  2.067  2.030  I 

1.047  radi— — I  1  1  1  1 

1  0.2  11.941  1.902  1.877  1.88  12.22  2.175  2.140  2.12  12.038  2.38  2.38  2.32  13.010  2.98  2.913  2.880  I 

! - 1  111: 

1  0.8  12.480  2.430  2.404  2.38  12.815  2.78  2.78  2.721  13.378  3.310  3.200  3.233  13.803  3.7M  3.734  3.0%  I 

1  0.2  10.819  0.808  0.792  0.784  10.930  0.917  0.98  0.8%  11.113  1.090  1.070  1.005  11.272  1.2(7  1.28  1.218  I 

i - 1  1  I  1  I 

90.0  0.40  11.487  1.437  1.437  1.423  11.700  1.005  1.043  1.027  12.02  1.979  1.98  1.98  12.80  2.203  2.28  2.80  1 

1.38  radi - 1  till 

1  0.8  12.107  2.003  2.037  2.017  11409  181  1321  1300  11803  126  1708  1740  13.84  3.207  3.104  3.18  I 

i - 1  111! 

!  0.2  11099  1044  1009  ISD  13.0%  3.023  198  198  13.007  3.23  3.545  3.309  14.193  4.12  4.08  4.012  1 


I  0.20  10.231  0.372  1.372  1.312  !1.182  1.901  1.9n  1.M5  il.SB2  2.399  2.09  2.202  11.981  3.177  3.309  2.712  I 

- 1  I  I  !  i 

M.0  dqi  0.90  10.928  0.087  2.912  2.930  !2.189  3.327  3.951  3.097  !2.93D  9.722  9.088  9.079  i3.979  5.921  6.129  5.115  ! 

1.097  r^i - 1  III! 

I  0.90  10.913  0.989  9.173  3.985  13.139  3.0SB  5.239  9.370  19.202  9.7n  7.009  5.850  13.299  8.990  8.789  7.3S  ! 

I - 1  I  I  I  I 

I  0.80  10.792  1.277  3.392  9.300  19.009  9.SS2  9.792  3.993  15.927  8.799  9.002  7.SBS  19.806  10.99  11.33  9.973  ! 


I  0.20  10.289  0.993  1.997  1.991  11.978  L3B2  2.999  2.008  11.979  3.189  3.302  2.735  12.982  3.999  9.190  3.9B  1 

I - 1  I  I  I  I 

90.0  *|l  0.90  10.335  0.892  3.993  3.090  12.738  9.912  9.511  3.812  13.999  5.907  9.115  5.103  19.397  7.907  7.998  9.399  I 

1.571  ridl - 1  1  1  I  I 

1  0.90  10.797  1.Z37  5.229  9.3»  13.927  9.327  9.S5D  5.997  15.257  8.9n  8.799  7.318  19.592  10.92  10.99  9.177  1 

1 - 1  111! 

1  0.80  10.991  1.597  9.799  5.930  15.0n  8.171  8.959  7.090  19.789  10.93  11.32  9.951  18.313  13.71  19.20  11.B  1 


C. 3.2(b)  Smedley-Wordsworth  Chord  SCF's  for  T-Joints 


The  Smedley-Wordsworth  chord  SCF's  for  T-joints  are  shown  on  the 
following  pages.  The  following  parameters  are  assumed  for  the 
Smedley-Wordsworth  figures: 

1)  -  =  2L/D  =  35.0 

The  Shell  d/D  limitations  have  not  been  imposed  for  the  SCF 
calculation. 


C-14 


WcriaarUi  Owiliy  SCF  Ciipitiftai 


T-joiit  Aiial  SCF  l>tM  Fkaitim 


6mm  «  110 

1  Gmm«  15.0 

Shm-  210 

1  1 

Ghm-  25.0 

!  0.25  !2.B(7  3.flSB  3.203  3.2BS  :2.02B  3.0B1  3.177  3.211  13.0M  3.148  3.180  3.106  13.207  3.228  3.206  3.145  t 

I - 1  I  I  i  * 

30.0  dql  0.S0  15.223  5.799  6.084  6.094  15.296  5.719  5.907  5.871  !5.492  5.734  5.793  5.678  15.716  5.823  5.775  1583  I 

0.524  radi - 1  it!! 

!  0.75  !1185  9.361  9.793  9.527  !1140  9.012  9.291  9.015  !1276  8.808  1891  1SS2  !1516  1807  1730  1308  t 

! - 1  II!! 

!  1.00  111.80  13.88  14.49  13.71  111.53  13.08  1146  12.75  !11.48  12.46  12.57  11.87  111.66  12.26  12.15  11.9  I 


!  0.25  !1194  1607  1961  4.29  !127B  1627  192  4.12  1142  1694  1920  4.12  11567  1774  1942  4.074  I 

! - 1  !  I  !  I 

410  dql  0.50  11771  6.724  7.433  7.909  1182  6.670  7.267  7.672  !6.106  6.714  7.165  7.466  16.354  6.821  7.157  7.364  I 

0.72  radi - 1  I  I  !  ! 

!0.2!172  10.47  11.54  12.2  !1829  10.21  11.11  11.54  19.059  10.11  10.2  11.10  19.39  10.16  10.67  10.87  1 

I - 1  II!! 

I  1.00  112.2  14.96  16.43  16.2  11119  14.36  15.9  15.9  !1133  13.96  14.87  1109  I116Z  119  14.S  14.66  I 


I  0.2  11218  1814  4.364  4.89  11294  1827  4.320  4.772  1149  18M  4.304  4.690  11557  199  4.318  4.651  I 


60.0  dql  0.9  11749 

1.047  rqi - 1 

!  0.2  18.52 


III! 

7.026  8.19  9.09  11869  7.001  7.99  18»  16.104  7.063  7.905  1631  16.345  7.12  7.903  1531  I 
III! 
10.60  12.27  1152  11643  10.47  11.2  1110  11933  10.46  11.2  12.2  !9.29  10.57  11.66  12.54  1 


!  1.00  111.57  14.63  16.9  112  111.65  14.30  16.29  17.65  111.94  14.14  llO  17.05  I12.S  14.19  1165  16.2  I 


!  0.2  13.084  3.89  4.609  5.39  13.19  3.862  4.547  112  11240  1892  4.511  1097  11343  121  4.59  102  I 

I - 1  I  !  !  ! 


2.0  dql  0.9  11569 

1.571  rqi - ! 

I  0.2  18.220 


7.205  8.62  10.02  1169  7.19  8.59  9.72  11842  7.19  1422  9.59  16.02  7.262  1397  9.443  I 


0.2  11220  10.81  1105  14.93  11315  10.67  12.2  14.51  11350  10.64  12.9  14.13  11814  10.71  12.41  1193  I 

- 1  III! 

1.2  Ill.M  14.76  17.76  20.10  111.14  14.45  17.22  19.45  111.39  14.9  16.2  1186  111.2  14.33  16.2  1154  I 


ArdSHorth-Siitf  cv  SCF  Co^Mtitin 

T-joiot  Aiial  ST  Saddle  PtBition 


12.0 


•  1S.0 


20.0 


210 


TheU 


Tan  >i  Beta  «d/0  BeU  «d/D  1  BeU  >d/D  Beta  -020 

t/T  !  0.3  i  0.5  I  0.7  !  0.9  !  0.3  1  0.5  !  0.7  !  0.9  1  0.3  i  0.5  I  0.7  I  0.9  !  0.3  !  0.5  !  0.7  5  0.9 


!  0.25  I0.S92 

0.973 

0.770 

! _ 

0.402  11.115 

1.216 

.  t 

0.963 

0.502  11.487 

— 

1.622 

1.2B4 

0.670  11.858 

2.03 

1.606 

0.808  i 

1 

1  1 

30.0  deoi  0.50  il.TM 

1.946 

1.541 

0.804  12.230 

2.433 

1.927 

1 

l.OOS  12.974 

1244 

2.569 

1.340  11717 

4.056 

1212 

1.676  ! 

0.5M  radi - ! 

1  0.75  !2.670 

»  -  » 

2.920 

2.312 

1.206  13.346 

f 

3.650 

2.890 

1.506  14.461 

4.867 

1654 

2.011  11576 

6.064 

4.818 

1 

1 

2.514  i 

1 

1  1 

;  1.00  13.569 

3.893 

3.083 

1.609  14.401 

4.867 

3.354 

2.011  11948 

6.489 

1139 

2.661  17.435 

8.112 

6.424 

1352  ! 

;  0.25  :i.6i8 

i.aoe 

i.SlO 

0.865  12.023 

2.260 

1.887 

1.081  12.698 

1014 

2.516 

1.442  11372 

1 

1768 

1146 

1.802  : 

■  J 

410  deal  0.50  13.237 

3.617 

3.020 

i 

1.730  14.047 

4.521 

3.775 

1 

2.163  11396 

6.028 

1033 

2.884  16.745 

7.536 

6.292 

1606  ! 

0.785  rad) - ! 

1  0.75  14.856 

1426 

4.530 

1 

1 

2.596  16.070 

f 

6.782 

1662 

1 

1 

1245  18.094 

1 

9.043 

7.550 

1 

4.326  110.11 

1 

11.30 

9.438 

5.406  1 

1  1.00  16.475 

7.234 

6.040 

1 

3.461  18.094 

9.043 

7.550 

4.326  110.79 

12.06 

10.06 

1768  11149 

15.07 

12.3 

7.211  1 

1  0.25  12.293 

2.598 

2.237 

1.354  12.867 

3.248 

2.797 

1.693  118Z3 

1 

4.331 

1729 

2.257  14.778 

t 

5.413 

4.662 

2.821  ! 

60.0  deoi  0.50  14.587 

1197 

4.475 

1 

2.709  11734 

6.496 

IW 

1 

1386  17.646 

8.662 

7.459 

4.515  19.557 

10.82 

9.324 

1643  ! 

1.047  radl - ! 

!  0.75  16.881 

7.795 

6.713 

1 

1 

4.063  18.601 

9.744 

8.391 

1 

1079  i  11.46 

12.99 

11.18 

1 

t 

6.772  114.33 
> 

16.24 

1198 

8.465  : 

1  i 

1.00  19,175 

10,39 

8.951 

5.418  ill. 46 

12,99 

11.18 

1 

6.772  11129 

17.32 

14.91 

9.030  119.11 

21.65 

18.64 

11.3  I 

1  0.25  12.937 

3.360 

2.958 

l.s6l  ‘3.671 

4.200 

3.697 

2.326  14.695 

1 

1600 

4.930 

1102  16.119 

1 

7.001 

6.162 

1678 

90.0  (tool  0.50  15.874 

6.721 

1916 

1 

J.723  17.343 

5.401 

/  •  v95 

1 

4.«3  19.790 

11.20 

9.860 

6.206  112.23 

14.00 

12.32 

7.756 

1.571  rad! - 1 

1  0.75  18.811 

10.08 

8. 874 

i 

5.584  111.01 

12.60 

11.09 

1 

6.980  114.66 

16.80 

14.79 

1 

1 

9.307  118.35 

1 

21.00 

18.48 

11.63 

t  i 

1  1.00  111.74 

13.44 

11.83 

7,446  ii4.68 

16. 80 

14.79 

1 

9.307  119.58 

22.40 

19.72 

12.41  124.47 

3.00 

24.65 

1151 

QOn  AVAILABLE  TO  DUC  DOBS  NOT  FESMIP  FOLLT  LEGIBLE  BEPRODUCTIO*, 


Mordao  th  awll»y  SCF  CafiiUtiOB 


T-jdnt  In-Plam  31  CroMi  Pontion 


6im>  12.0  6mh>  15.0 

Gaat*  20.0  6aai«  2L0 

Tau  ■!  Brta  «d/0  Biti  >d/0 

t/T  1  0.3  1  0.5  1  0.7  1  0.9  1  0.3  1  0.5  1  0.7  1  0.9 

BMa>d/D  Brti>d/D 

0.3  1  0.5  1  0.7  1  0.9  1  0.3  1  0.5  1  0.7  1  0.9 

30.0  dq 
0.524  rad 

0.25  10.607  0.791  0.946  1.044  10.694  0.906  1.081  1.194  10.825  1.075  1.285  1.419  10.943  1.229  1.469  1.623 

0.50  11.057  1.378  1.647  1.819  11.209  1.575  1.883  2.080  11.437  1.872  2.23B  2.471  11.612  2.141  2.558  2.826 

-  1  1  1  j 

0.75  11.462  1.906  2.278  2.516  11.672  2.179  2.604  2.877  11.987  3.095  3.419  12.272  1961  3.539  3.908 

'  '  1  1  1  i 

1.00  11.841  1399  1868  3.167  11105  1743  3.278  3.621  11501  3.260  3.896  4.303  11860  3.727  4.455  4.920 

:  0.25  10.865 

*  »  « 

1.009 

1.079 

1.066  10.989 

1 

1.153 

1.233 

1.219  11.175 

1.370 

1.466 

1.449  11.343 

1 

1.567 

1.676 

1.657  1 

45.0  dql  0.50  11.506 
0.785  rad! - 1 

1.756 

1.879 

1 

1.857  11.721 

1 

t 

1008 

1148 

1123  11046 

1 

1 

1386 

1553 

1 

1523  11339 

1 

I 

1728 

1918 

1885  1 

1  1  0.75  11083 

* 

1429 

1599 

1569  11381 

t 

1778 

1971 

1937  11830 

1 

3.301 

3.531 

3.490  13.235 

1 

3.774 

4.037 

3.991  1 

1 

1  1  1.00  11622 

3.058 

3.271 

1 

3.234  11998 

3.497 

3.740 

1 

3.697  13.562 

4.155 

4.445 

1 

4.m  14.073 

4.751 

5.082 

5.023  I 

1  0.25  11.063 

1  i  _ 1 

1.162 

1.165 

1.080  11.216 

1 

1.329 

1.332 

1.234  11.445 

1 

1.579 

1.583 

1.467  11.652 

I 

1.806 

1.810 

1.677  1 

1  1 

!  60.0  dql  0.50  11.852 

1  1.047  radi - 1 

1024 

1029 

t 

1.880  11117 

1 

1 

1314 

1320 

1 

1149  11516 

1 

1 

1750 

1757 

1 

1554  11876 

1 

3.144 

3.152 

1920  1 

1  1  0.75  11561 

» - 1 

1800 

1807 

1600  11928 

1 

3.201 

3.209 

1973  13.480 

1 

3.804 

3.314 

3.533  13.979 

1 

4.350 

4.360 

4.039  1 

1  '■■■■  ■  1 

1  1.00  13.224 

3.525 

3.533 

( 

3.273  13.686 

4.030 

4.040 

1 

3.742  14.381 

4.789 

4.801 

1 

4.448  15.009 

5.475 

5.489 

5.085  1 

!  0.25  !1.231  1.286  1.231  1.089  !1.40e  1.470  1.407  1.245  11.673  1.747  1.672  1.480  11.913  1.997  1.912  1.692 

»  - I  f  I  ( 

I  J  I  I  I 

90.0  dcq:  0.50  12.144  2.239  2.143  1.896  12.452  2.559  2.450  2.168  12.914  3.042  2.912  2.576  13.331  3.478  3.29  2.946 

1.571  ridl - 1  1  1  1 


Uordaortli-Sitdlfy  SCF  CoipiiUtiaa 


T-joint  Out-of-flint  SCF  Saddlo  Position 


'  i  1 

•  I  I 

SlHI  I 

12.0 

1 

1 

SaMt* 

15.0 

4 

1 

Baas  « 

20.0 

t 

1 

5mm  * 

25.0 

!  TliiU «  !  Tai  «|  Brti  -d/D  !  Brtt  >d/D  Brta  *d/D 

■  1  t/T  !  0.3  {  0.5  0.7  I  0.9  !  0.3  I  0.5  !  0.7  !  0.9  1  0.3  i  0.5  !  0.7  1  0.9  !  0.3 

BlU  ad/D 

1  0.5  1  0.7 

1  0.9 

!  0.2S  !0.529 

0.773 

0.825 

0.556  10.662 

1 

r 

0.967 

t 

1.031 

/  / 

0.695  !0.883 

1 

1.289 

1.375 

0.927  11.103 

1.612 

1.719 

1.159 

;  30.0  dqi  0.50  il.0S9 

0.529  rid! - 1 

1.547 

1.650 

1 

1.112  11.324 

1 

1 

1.934 

2.062 

1 

1.390  !1.766 

» 

1 

2.579 

2.750 

1 

1.854  12.207 

! 

3.224 

3.438 

2.318 

!  0.75  I1.SB9 

2.321 

2.475 

1.669  !1.986 

1 

2.902 

3.094 

2.066  !2.649 

3.869 

4.125 

2.781  13.311 

4.837 

5.157 

3.477 

;  1.00  12.119 

3.095 

3.300 

1 

2.225  12.649 

3.869 

4.125 

1 

2.781  !3.532 

5.159 

5.500 

1 

3.709  14.415 

6.449 

6.876 

4.636 

!  0.25  !0.S72 

1.347 

1.561 

1.176  11.090 

1 

1.684 

1.951 

1.470  !1.4S4 

2.245 

2.602 

1.960  11.818 

2.807 

3.252 

1450  1 

45.0  dqi  0.50  .'1.745 
0.785  ridl - -1 

2.694 

3.122 

1 

2.352  12.181 

1 

1 

3.368 

3.903 

1 

2.940  12.906 

1 

1 

4.491 

5.204 

1 

3.920  13.636 

I 

1 

5.614 

6.506 

4.900  1 

!  !  0.75  !2.61B 

4,042 

4.683 

3.528  !3.272 

1 

5.053 

5.854 

4.410  14.363 

6n737 

7.806 

5.880  15.454 

8.421 

9.757 

7.351  1 

!  1.00  13.490 

5.389 

6.245 

1 

4.704  14.363 

6.737 

7.806 

1 

5.880  15.817 

8.983 

10.40 

1 

7.841  17.272 

11.22 

13.01 

9.801  ; 

i  !  0.25  11.168 

1.863 

2.267 

1.822  11.460 

1 

2.329 

2.833 

2.278  !1.947 

3.106 

3.778 

3.037  12.434 

3.8BZ 

4.7Z3 

3.796  1 

60.0  dcgl  0.50  12.337 

1  1.047  rad! - ! 

3.727 

4.534 

1 

3.644  12.921 

1 

1 

4.659 

5.667 

1 

4.556  !3.B95 

1 

6.212 

7.557 

$ 

6.074  14.868 

1 

1 

7.765 

9.446 

7.593  ! 

I  0.75  !3.5Q6 

r  f - _  « 

5.591 

6.801 

5.467  !4.382 

f 

6.989 

8.501 

6.834  15.842 

9.318 

11.33 

9.112  17.303 

11.64 

14.16 

11.39  1 

;  1.00  !4.674 

7.455 

9.068 

7.219  !S.842 

9.318 

11.33 

1 

9.112  !7.790 

12.42 

15.11 

r 

12.14  19.737 

15.53 

18.89 

15.18  1 

:  0.25  !1.437 
.  »  1  * 

2.346 

2.954 

2.486  11.796 

1 

2.932 

3.692 

3.108  12.395 

3.910 

4.923 

4.144  12.994 

4.887 

6.154 

5.180  ! 

90.0  dql  0.50  i2.874 
1.571  rad! — — ! 

4.692 

5.908 

1 

4.973  13.593 

* 

1 

5.865 

7.385 

1 

6.216  !4.791 

1 

1 

7.820 

9.847 

t 

8.288  15.989 

1 

1 

9.775 

12.30 

10.36  1 

1  0.75  J4.312 

J  —  < 

7.038 

8,862 

7.459  !5.390 

1 

8.797 

11.07 

9.324  17.187 

11.73 

14.77 

12.43  18.984 

14.66 

18.46 

15.54  1 

:  1.00  15.749 

9.384 

11.81 

9.946  17.187 

11.73 

14.77 

12.43  19.583 

15.64 

19.69 

16.57  111.97 

19.55 

24.61 

20.72  ! 

K-Joilit  Md  SEF 
OMSMiof  Mi 


6Mh>  12.0  1 

1  SHf  ao 

8m-  20lO  1 

1  1 

1  Gm-  25.0 

I  S.2D  :0.4Z7 
I - 1 

30.0  *|>  0.40  10.918 

0.321  radi - 1 

I  0.00  !1.<37 


0.414  0.400  0.400  iO.495  0.4n  0.4n  0.404  10.000  0.3B  0.371  0.302  :0.0i0  0.070  0.002  0.1SS  I 

.'111 

0.891  0.873  0.800  11.005  1.034  1.013  0.991  11.290  1.232  1.227  1.209  11.497  1.453  1.424  1.403  ! 

I  I  I  I 

1.394  1.307  1.347  11.007  1.018  1.SB0  1.302  S2.019  1.931  1.921  1.892  12.313  2.273  2.229  2.190  I 


I - 1  I  i  i  • 

I  0.80  11.974  I.91S  1.878  1.830  11290  1222  1179  1147  11774  1092  1039  lOOO  13.219  3.123  3.002  3.017  ! 


!  0.20  10.724  0.702  0.0a  0.078  10.810  0.813  0.799  0.787  11.017  0.987  0.907  0.953  11.180  1.143  1.122  1.100  1 

I - 1  Jill 


!  410  Mi  0*40  11S50 

1  0.783  rail - 1 

1  1  0.00  11431 

1  I - 1 

1  1 180  11345 


1  0.20  10.985  0.950  0.937  0.9Z3  11.143  1.109  1.087  1.0n  11.384  1.343  1.317  1.297  11.000  1.359  1.328  1.305  1 

1  1 - 1  i  I  I  i 

I  00.0  M>  0.40  11118  1005  1015  1.983  11487  1381  1337  1305  11970  1888  180  1789  13.453  3.351  3.285  3.230  1 
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! - 1  111! 
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Itef  97  DapaUtioi 

K-jaiit  b-PlMt  97 
Qiarri  Sidi  of  lUd 
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1  !  1 
3.087  3.13  13.273  3.375  3.444  3.494  13.543  3.473  1748  1805  I 
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1 
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1 
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5.20? 
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1 
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1.371  radi - 1  I  !  ' 
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NbrdMrtit-SiMUtr  SCF  Coipotatiin 
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I 
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0.811  11.229 

( 

1.356 

» 
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30.0  d^l  0.50  !1.525 
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1 
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1 
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1 
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OOP?  AVAILABLE  TO  OTIC  DOES  MOT  FEEMII  FULLY  LEGIBLE  REPBODUCTIOH 
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C.4  FINITE  ELEMENT  ANALYSES  RESULTS 


Finite  element  analyses  were  performed  on  the  connections  between  the 
column  tops  and  upper  hull  girders  and  between  the  corner  columns  and 
tubular  bracing  of  the  column-stabilized,  twin-hulled 
semi submersible.  The  overall  geometry  and  locations  of  the  two 
connections  are  indicated  in  Figures  C.4-1  and  C.4-2.  Longitudinal 
and  transverse  girders  (8.2  m  deep)  coincide  with  the  column  faces. 
The  columns  are  10.6  by  10.6  m  in  cross  section. 

Some  dimensions  of  interest  are: 

Overall  length  96.0  m 

Overall  width  65.0  m 

Lower  Hulls  (two) 

Length 
Width 
Depth 

Stability  Columns  (six) 

Size  (square  w/  rounded  columns)  10.6x10.6  m 
Transverse  spacing  (center- to-center)  54.0  m 
Longitudinal  spacing  33.0  m 

Upper  Hull 

Length  77.0  m 

Width  65.0  m 

Depth  8.2  m 

C.4.1  Column-Girder  Connection 


96.0  m 
16.5  m 
8.0  m 


The  location  of  the  connection  is  shown  in  Figure  C.4-1.  The  joint 
dimensions  are  given  in  Figure  C.4-3.  The  loading  analyzed  was  a 
combined  axial,  shear  and  moment  load.  The  SCF  is  defined  as: 
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MAXIMUM  PRINCIPAL  STRESS 


SCF . . 

NOMINAL  STRESS  IN  GIRDER  (P/A  +  M/S) 

The  moment  M  is  due  to  a  combination  of  moment  and  shear  load. 

The  maximum  SCF  was  found  In  the  gusset  plate  connecting  the 
transverse  girder  and  column  top,  at  the  edge  of  the  gusset  plate  In 
the  weld  between  the  gusset  web  and  flange.  It  was  equal  to  1.66.  The 
SCF  1n  the  longitudinal  girder  at  the  windlass  cutouts  reached  a 
value  of  1.87.  Figure  C.4-4  shows  the  equivalent  stress  variation 
over  the  entire  connection.  The  maximum  stress,  as  already  noted, 
occurs  In  the  crotch  region.  Figure  C.4-5  shows  an  equivalent  stress 
contour  plot  of  the  windlass  holes.  Table  C.4-1  summarizes  the  SCFs. 
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MEMBER  LOCATION  DIRECTION  SCF 

TO  WELO 


Center  Column  Middle  of  Parallel  1.66 

Transverse  Gusset 

Girder-Column  Gusset-Girder  Perpendicular  1.37 

Connection  Connection 

2. 3x2. 3x1.1  m  Gusset-Colunm  Perpendicular  1.10 

Gusset  Connection 


Girder  Flange  Parallel  1.05 


Longitudinal  Middle  of  Parallel  1.52 

Glrder-Colunm  Gusset 


Connection  Gusset-Girder  Perpendicular  1.14 

Connection 

1.1x1. lx. 55  m  Gusset-Colunm  Perpendicular  1.05 

Connection 


Girder  Flange  Parallel  1.01 


Exterior  Bottom  Right  Parallel  1.87 

Longitudinal  Corner  of 

Girder  Exterior  Hole 


9  Windlass  Upper  Left  Parallel  1.73 

Holes  Corner  of 

Interior  Hole 


Table  C.4-1  Sinmary  of  SCFs  for  a  Coluan  ^  Girder  Connection 


Flgura  C.4-1  Overall  GaoMtry  of  vassal  and  Location  of 
Colusn-to-Girdar  connaetlon 
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Flours  c.4-2  Overall  GaoMtry  of  Vassal  and  Location  of 
Cornar  coluan  Tubular  connaction 
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Pigur.  c.4-3  Dimmnmionm  of 

coluan-to-Glrdor  Coanoction 


Figure  C.4-4  Equivalent  Stress  Contour  Plot  for  the 
Coluan>to-Girder  Connection 
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-5  Equivaiant  Strass  Contour  Plot  of  Windlass  Rolss 
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CL 
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D 

*^tot 

h 

^2 

E 

H 

I 

I 

lo 

K 

Ks 

L 

N 

Re 

S 

S 

SCF 

St 

Si 

T 


'^max 

'^mln 


’m 

Yr 


Coefficient  of  drag 

Design  lift  coefficient 

Base  lift  coefficient 

Fatigue  damage 

Total  fatigue  damage 

Fatigue  damage  due  to  vortex  shedding 

Fatigue  damage  due  to  storm 

Modulus  of  elasticity 

Submerged  length  of  member 

Member  moment  of  Inertia 

Turbulence  parameter 

Turbulence  parameter 

Constant  representing  member  fixity 

Stability  parameter 

Span  between  member  supports 

Number  of  cycles  to  failure  at  hot  spot  stress  range 

Reynolds  number 

Hot  spot  stress  range 

Member  section  modulus 

Stress  concentration  factor 

Strouhal  number 

Corresponding  hot  spot  stress  range 
Wave  period 

Time  for  which  '^mln  exceeded 

Flow  velocity  normal  to  member  axis 

Maximum  orbital  velocity  due  to  wave  motion 

Maximum  water  particle  velocity 

Minimum  required  for  motion 

Reduced  velocity 

Member  midspan  deflection 

Maximum  member  midspan  deflection 

Refined  maximum  member  midspan  deflection 


**1 

y(x) 

y'(x) 


e 

V 


’’N 

P 


Maximum  modal  ^lltude 
Natural  frequency  coefficient 
Pit  depth 
Member  diameter 

Member  vortex  shedding  frequency 

Turbulence  parameter 

Member  bending  stress 

Member  maximum  bending  stress 

Member  maximum  hot  spot  stress 

Member  natural  frequency 

Member  vortex  shedding  frequency 

Mass  of  member  per  unit  length  excluding  marine  growth 

Effective  mass  per  unit  length 

Mass  of  meirt>er  per  unit  length  Including  marine  growth 

Generalized  mass  per  unit  length  for  mode  J 

Mode  of  vibration 

Member  end  condition  coefficient 

Total  number  of  occurrences  per  year 

Actual  number  of  cycles  at  hot  spot  stress  range 

Number  of  oscillations  during  one  wave  cycle 

Nominal  caisson  thickness 

Applied  velocity 

Critical  wind  velocity 

Reduced  velocity 

Load  per  unit  length 

Weight  per  unit  length  of  member 

Weight  per  unit  length  of  supported  Items 

Fundamental  mode  shape 

Equivalent  fundamental  mode  shape 

Logarithmic  decrement  of  damping 

Damping  ratio 

Kinematic  viscosity 

Ratio  of  midspan  deflection  to  member  diameter  (Y/d) 
Mass  density  of  fluid 
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VORTEX  SHEDDING 


D.l.  INTRODUCTION 


When  a  fluid  flows  about  a  stationary  cylinder,  the  flow  separates, 
vortices  are  shed,  and  a  periodic  wake  Is  foimed.  Each  time  a  vortex 
Is  shed  from  the  cylinder,  the  local  pressure  distribution  Is 
altered,  and  the  cylinder  experiences  a  time-varying  force  at  the 
frequency  of  vortex  shedding. 

In  steady  flows,  vortices  are  shed  alternately  from  either  side  of 
the  cylinder  producing  an  oscillating  lift  force  transverse  to  the 
flow  direction  at  a  frequency  equal  to  that  at  which  pairs  of 
vortices  are  shed.  In  the  flow  direction.  In  addition  to  the  steady 
drag  force,  there  Is  a  small  fluctuating  drag  force  associated  with 
the  shedding  of  Individual  vortices  at  a  frequency  twice  that  of  the 
lift  force. 

As  the  flow  velocity  Increases,  the  vortex  shedding  frequency 
Increases.  Thus,  provided  the  flow  velocity  Is  high  enough,  a 
condition  will  be  reached  where  the  vortex  shedding  frequency 
coincides  with  the  natural  frequency  of  the  flexible  element. 

In  general,  marine  members  and  appurtenant  pipework  are  of  a  diameter 
and  length  that  preclude  the  occurrence  of  In-line  vibrations  Induced 
by  vortex  shedding.  However,  all  susceptible  members  must  be 
analyzed  to  ensure  that  the  stresses  due  to  In-line  vibrations  and 
possible  synchronized  oscillations  are  small  and  do  not  result  In  a 
fatigue  failure. 

Response  to  vortex  shedding  cannot  be  predicted  using  conventional 
dynamic  analysis  techniques  since  the  problem  Is  non-linear.  The 
motion  of  the  structure  affects  the  strength  of  the  shedding  which. 
In  turn  affects  the  motion  of  the  structure.  This  feedback  mechanism 
causes  the  response  to  be  either  significantly  large  or  negligibly 
small.  Once  excited,  there  Is  also  a  tendency  for  the  vortex 
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shedding  frequency  to  synchronize  with  the  natural  frequency  of  the 
structure.  This  results  In  sustained  resonant  vibration  even  If  the 
flow  velocity  moves  away  from  the  critical  velocity. 

Oscillations  can  be  predominantly  In-line  with  the  flow  direction  or 
transverse  to  It.  In-line  motion  occurs  at  lower  flow  velocities 
than  transverse  or  cross-flow  motion,  but  the  latter  Is  Invariably 
more  severe  and  can  lead  to  catastrophic  failure  due  to  a  small 
number  of  cycles  of  oscillation. 

Response  to  vortex  shedding  Is  further  complicated  as  the 
excitatlonal  force  Is  not  necessarily  uniform  along  the  length  of  the 
members  and  the  actual  amplitude  of  oscillation  depends  to  a  large 
extent  on  the  degree  of  structural  damping. 

D.2.  VORTEX  SHEDDING  PARAMETERS 

A  number  of  parameters  are  common  to  this  phenomenon: 

Reduced  velocity  (V^) 

V  • 

where: 

V 

^n  = 

d  » 

Reynolds  number  (Rg) 

Rg  *  Vd/  V 

where: 

V  *  kinematic  viscosity  of  the  fluid 

The  Strouhal  number  (S^)  is  a  function  of  the  Reynolds  number  for 
circular  members.  The  Reynolds  number  for  typical  cylindrical 

members  under  storm  current  ranges  from  3.5  x  10®  to  1.0  x  10®.  The 
Strouhal  number  1s  reasonably  approximated  as  0.21  for  this  range  of 
Reynolds  numbers. 


V/f„d 

flow  velocity  normal  to  the  member  axis 
fundamental  frequency  of  the  member  (H^) 
diameter  of  the  member 
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Vortex  Shedding  Frequency  (f^) 


=  — g - *  vortex  shedding  frequency  of  the  member 

If  the  vortex  shedding  frequency  of  the  member  coincides  with  the 
natural  frequency  of  the  member,  resonance  will  occur. 


Stability  parameter  (K^) 


•^S 

s 

2 

2  m  «  /  p  d 

where: 

6 

= 

2  ir  e  =  logarithmic  decrement 

e 

s 

damping  ratio 

P 

a 

mass  density  of  the  fluid 

m 

3 

effective  mass  per  unit  length 

s 

/o(ni)ly(x)l^dx 

Jo  iy'(«)i^<ix 

L 

3 

span  between  member  supports 

m 

- 

mass  of  member  per  unit  length 

y(x). 

&  y(x) 

fundamental  mode  shapes  as  a  function  of  the 
ordinate  x  measured  from  the  lower  support 
along  the  longitudinal  axis  of  the  mend>er 

As  given  In  References  0.1  and  0.2,  the  effective  mass  Is  used  to 
equate  the  real  structure  with  an  equivalent  structure  for  which 
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deflection  and  stability  parameters  are  known.  The  deflected  form  of 
this  equivalent  structure  Is  a  cantilever,  while  typical  structure 
members  and  appurtenances  deflect  as  a  simply  supported  beam.  Hence, 
the  equivalent  structure  has  a  mode  shape  given  by: 

y'(x)  =  a  -  a  cos(^) 

while  the  real  structure  has  a  mode  shape  given  by: 
y(x)  =  a  s1n(^) 

Substituting  Into  the  effective  mass  formulation,  we  obtain: 

/glmlfa  sin  -pl^dx 
/q  (a  -  a  cos 


where: 


a  »  maximum  modal  amplitude 

Integration  of  the  above  equation  leads  to  the  relationship: 
m  =  2.205  m  for  simple  supported  span 

m  »  1.654  m  for  fixed  supports 

m  »  m  for  cantilever  span 

Dampinq  Ratio 

Welded  marine  structures  exhibit  very  low  values  of  structural 
damping.  Vibratory  energy  Is  typically  dissipated  by  material  and 
aerodynamic  (radiation)  damping.  Individual  members  subjected  to 
large  vibratory  motions  dissipate  energy  through  the  connections  to 
the  main  structure  largely  as  dispersive  bending  and  compression 
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waves.  When  only  Isolated  members  undergo  large  vibration  response, 
energy  dispersion  exceeds  reflected  energy  and  represents  a  major 
source  of  damping. 

Structural  members  may  be  grouped  Into  two  classes,  depending  on  the 
fixity  of  their  supports.  Tubular  braces  welded  on  to  regions  of 
high  rigidity,  such  as  structure  columns  or  legs,  are  defined  as 
Class  1  members.  Tubular  braces  welded  on  to  regions  of  low 
rigidity,  such  as  other  braces,  are  defined  as  Class  2  members.  The 
damping  ratio  applicable  for  structural  members  are: 

Structural  Member  -  Class  1  Damping  ratio  e  «  0.0035 
Structural  Member  -  Class  2  Damping  ratio  e  =  0.0015 

Although  the  recommended  damping  ratios  are  for  vibrations  In  air, 
they  may  be  conservatively  used  for  vibrations  In  water. 

Non-structural  continuous  members,  such  as  tubulars  supported  by 
multiple  guides,  have  both  structural  and  hydrodynamic  damping.  The 
hydrodynamic  damping  occurs  due  to  sympathetic  vibration  of  spans 
adjacent  to  the  span  being  evaluated  for  shedding.  Recent  work  by 
Vandiver  and  Chung  (Reference  0.3)  supports  the  effectiveness  of 
hydrodynamic  damping  mechanism.  The  lower  bound  structural  damping 
ratio  for  continuous  tubulars  supported  by  loose  guides  Is  given  as 
0.009  by  Blevins  (Reference  0.4).  The  applicable  damping  ratios  are 
assumed  to  be: 

Non>Structural  Members  -  Continuous  Spans 

Damping  ratio  c  -  0.009  In  air 
Damping  ratio  e  >  0.02  In  water 

Natural  Frequency 

The  fundamental  natural  frequency  (In  Hz)  for  uniform  beams  may  be 
calculated  from: 
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^n 

s 

^  (EI/BjlV 

where: 

I 

s 

the  moment  of  Inertia  of  the  beam 

^n 

« 

3.52  for  a  beam  with  fix-free  ends  (cantilever) 

a 

9.87  for  a  beam  with  pin-pin  ends 

a 

15.4  for  a  beam  with  fix-pin  ends 

a 

22.4  for  a  beam  with  fix-fix  ends 

L 

a 

length 

n 

= 

mode  of  vibration 

"1 

a 

mass  per  unit  length 

The  ainount  of  member  fixity  assumed  In  the  analysis  has  a  large 
effect  on  vortex  shedding  results,  because  of  Its  Impact  on  member 
stiffness,  natural  period,  amplitude  of  displacement,  and  member 
stress.  Hence,  careful  consideration  should  be  given  to  member  end 
conditions.  Members  framing  Into  relatively  stiff  members  can 
usually  be  assumed  to  be  fixed.  Other  members,  such  as  caissons  and 
risers,  may  act  as  pinned  members  If  supports  are  detailed  to  allow 
member  rotation. 

For  members  with  non-uniform  spans,  complex  support  arrangements  or 
non-uniform  mass  distribution,  the  natural  frequency  should  be 
determined  from  either  a  dynamic  analysis  or  from  Tables  provided  In 
References  0.5  and  0.6.  Reid  (Reference  0.7)  provides  a  discussion 
and  a  model  to  predict  the  response  of  variable  geometry  cylinders 
subjected  to  a  varying  flow  velocities. 

The  natural  frequency  of  a  member  Is  a  function  of  the  lumber's 
stiffness  and  mass.  For  the  purposes  of  vortex  shedding  analysis  and 
design,  the  member's  stiffness  properties  are  computed  from  the 
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member's  nominal  diameter  and  thickness.  The  member  mass  per  unit 
length  m  Is  taken  to  Include  the  mass  of  the  member  steel  Including 
sacrificial  corrosion  allowance,  anodes,  and  contained  fluid.  For 
the  submerged  portion  of  the  member,  the  added  mass  of  the 
surrounding  water  Is  also  Included.  This  added  mass  Is  the  mass  of 
water  that  would  be  displaced  by  a  closed  cylinder  with  a  diameter 
equal  to  the  nominal  member  outside  diameter  plus  two  times  the 
appropriate  marine  growth  thickness. 

Because  of  Insufficient  knowledge  of  the  effect  of  marine  growth  on 
vortex  shedding,  the  member  diameter  "d"  in  vortex-shedding 
parameter'i  V^,  Rg,  K^,  and  the  member  effective  mass  m  In  parameter 
Kj  do  not  Include  any  allowance  for  the  presence  of  marine  growth. 

D.3.  SUSCEPTIBILITY  TO  VORTEX  SHEDDING 

The  vortex  shedding  phenomena  may  occur  either  In  water  or  In  air. 
The  susceptibility  discussed  and  the  design  guidelines  presented  are 
applicable  for  steady  current  and  wind.  Wave  Induced  vortex 
shedding  has  not  been  Investigated  In  depth.  Since  the  water 
particle  velocities  In  waves  continually  change  both  In  magnitude  and 
direction  (I.e.  restricting  resonant  oscillation  build-up).  It  may  be 
reasonable  to  investigate  current- Induced  vortex  shedding  and 
overlook  wave  actions. 

To  determine  susceptibility  of  a  member  to  wind-  or  current- Induced 
vortex  shedding  vibrations,  the  reduced  velocity  (V^)  Is  computed 
first.  For  submerged  members,  the  stability  parameter  (K^)  is  also 
calculated.  Vortex  shedding  susceptibility  defined  here  Is  based 
upon  the  method  given  In  Reference  D.8,  with  a  modified  lower  bound 
for  current- Induced  shedding  to  reflect  present  thinking  on  this 
subject  (Reference  0.9). 
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D.3.1  In-Line  Vortex  Shedding 

In-line  vibrations  In  wind  and  current  environments  may  occur  when: 

Current  Environment  Wind  Environment 

1.2  s  <  3.5  1.7  <  V^<  3.2 

and  <1.8 

The  value  of  may  be  more  accurately  defined  for  low  values  from 
Figure  0-1,  which  gives  the  reduced  velocity  necessary  for  the  onset 
of  In-line  motion  as  a  function  of  combined  mass  and  damping 
parameter  (I.e.  stability  parameter).  Corresponding  amplitude  of 
motion  as  a  function  of  Is  given  on  Figure  0-2.  As  Illustrated  on 
this  Figure,  In-line  motion  Is  completely  supressed  for  values 
greater  than  1.8. 

Typical  marine  structure  members  (I.e.  braces  and  caissons  on  a 
platform)  generally  have  values  of  greater  than  1.8  In  air  but 
less  than  1.8  In  water.  Hence,  In-line  vibrations  with  significant 
amplitudes  are  often  likely  In  steady  current  but  unlikely  In  wind. 

0.3.2  Cross-Flow  Vortex  Shedding 

The  reduced  velocity  necessary  for  the  onset  of  cross-flow  vibrations 
In  either  air  or  In  water  Is  shown  on  Figure  0-3  as  a  function  of 
Reynold's  number,  Rg,  cross  flow  vibrations  In  water  and  In  air  may 
occur  when: 

Current  Environment  Wind  Environment 


3.9  s  <  9  4.7  <  <  8 

and  <  16 

The  cross-flow  v1brv*t<ons  of  members  in  steady  current  will 
Invariably  be  of  large  amplitude,  causing  failures  after  small  number 
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of  cycles.  Thus,  the  reduced  velocity  necessary  for  the  onset  of 
cross-flow  vibrations  In  steady  current  should  be  avoided. 

D.3.3  Critical  Flow  Velocities 


The  criteria  for  determining  the  critical  flow  velocities  for  the 
onset  of  VIV  can  be  expressed  In  terms  of  the  reduced  velocity 
(Section  D.2): 

Vcr  -  (Vr)cr  C'n-  '*) 
where: 

(Vr)  cr  *  In-line  oscillations  In  water 

»  1.7  for  In-line  oscillations  in  air 

*  3.9  for  cross-flow  oscillations  In  water 

*  4.7  for  cross-flow  oscillations  in  air 

0.4.  AMPLITUDES  OF  VIBRATION 

Amplitudes  of  vibrations  can  be  determined  by  several  methods.  A  DnV 
proposed  procedure  (Reference  0.8}  is  simple  to  apply  and  allows 
determination  of  member  natural  frequencies,  critical  velocities  and 
maximum  amplitudes  of  vortex-shedding  induced  oscillations.  The 
procedure  yields  consistent  results,  comparable  to  the  results 
obtained  by  other  methods,  except  for  osci Hatton  amplitudes.  The 
DnV  calculation  of  oscillation  amplitudes  Is  based  on  a  dynamic  load 
factor  of  a  resonant,  damped,  single-degree-of -freedom  system,  this 
approach  1s  not  valid  unless  the  nonlinear  relationship  between  the 
response  and  damping  ratio  Is  known  and  accounted  for.  Consequently, 
In-line  and  cross-flow  vortex  shedding  amplitudes  are  assessed 
separately. 
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D.4.1  In-Line  Vortex  Shedding  Amplitudes 


The  reduced  velocity  and  the  amplitude  of  vibrations  shown  on  Figures 
0-1  and  0-2,  respectively,  as  functions  of  stability  parameter  are 
based  on  experimental  data.  The  experimental  data  obtained  are  for 
the  cantilever  mode  of  deflection  for  in-line  and  cross-flow 
vibrations. 

Sarpkaya  (Reference  0.10)  carried  out  tests  on  both  oscillatory  flow 
and  uniform  flow  and  observed  smaller  amplitudes  of  vibration  for  the 
oscillatory  flow  than  for  the  uniform  flow.  It  is  also  suggested  by 
King  (Reference  0.1)  that  the  maximum  amplitude  for  an  oscillatory 
flow  is  likely  to  occur  at  a  value  in  excess  of  1.5  (as  opposed  to 
1.0  assumed  by  OnV)  and  that  an  oscillation  build-up  of  about  15 
cycles  is  required  before  "lock-in”  maximum-amplitude  vibration 
occurs.  In  light  of  this  evidence,  the  amplitude  of  vibrations  shown 
in  Figure  0-2  is  based  on  Hal  1am  et  a1  (Reference  0.2)  rather  than 
the  OnV  (Reference  0.8). 

Since  typical  marine  structure  members  have  stability  parameters  (K^) 
in  excess  of  1.8,  in-line  vibrations  of  these  members  in  air  are 
unlikely. 

0.4.2  Cross-flow  Vortex  Shedding  Amplitudes 

The  amplitude  of  the  induced  vibrations  that  accompanies  cross-flow 
vibration  are  generally  large  and  creates  very  high  stresses. 
Therefore,  it  is  desirable  to  preclude  cross-flow  induced 
vibrations.  Figure  0-4  illustrates  a  curve  defining  the  amplitude  of 
response  for  cross-flow  vibrations  due  to  current  flow  and  based  on  a 
cantilever  mode  of  deflection. 

Cross-flow  oscillations  in  air  may  not  be  always  avoidable,  requiring 
the  members  to  have  sufficient  resistance.  The  OnV  procedure 
(Reference  0.8)  to  determine  the  oscillation  amplitudes  is  derived 
from  a  simplified  approach  applicable  to  vortex  shedding  due  to 
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steady  current,  by  substituting  the  mass  density  of  air  for  the  mass 
density  of  water.  Hence,  the  oscillation  amplitude  Is  not  linked 
with  the  velocity  that  causes  vortex- Induced  motion.  The  resulting 
predicted  amplitudes  are  substantially  higher  than  amplitudes 
predicted  based  on  an  ESOU  (Reference  D.ll)  procedure  that  accounts 
for  Interaction  between  vortices  shed  and  the  forces  Induced. 

The  Iterative  ESOU  procedure  to  determine  the  amplitudes  can  be 
simplified  by  approximating  selected  variables.  The  peak  amplitude 
Is  represented  In  Equation  9  of  the  ESOU  report  by: 

Y  0.00633  (?  1 

I  “  ■  ""1 - ^  mT  ^ 

J  ^t 

Using  this  formulation,  a  corresponding  equation  can  be  established 
for  a  structure,  while  making  assumptions  about  the  Individual 
parameters.  Following  step  3  of  the  procedure,  the  parameters  may  be 
set  as: 


0.0795  CL. 
- 


0-11 


base  lift  coefficient 

3 

0.29  high  Reynolds  number 

9 

0.42  low  Reynolds  number 

CLj  - 

design  lift  coefficient 

* 

CLo  X  f„  X  lo  x-jj  X  1.2 

^ar  “ 

turbulence  parameter 

9 

1.0 

I 

turbulence  parameter 

- 

1.0 

'o  - 

turbulence  parameter 

. 

0.45 

Evaluating  the  equation  based  on  the  high  Reynolds  number  (Re  > 
500,000)  leads  to: 


»  0.0795  (0.29)(1.0)(0.45)(1.0)(1.2)/[Kj(0.2)^l 
0.3114 


'’N 


(high  Reynolds  number.  Re  >  500,000) 


or 


»  Oi^SlO  Reynolds  number.  Re  <  500,000) 


The  amplitude  can  also  be  determined  Iteratively  by  utilizing  the 
ESDU  recommended  turbulence  parameter  and  following  steps  1  through 
5. 

Step  1:  Determine  correlation  length  factor,  1^.  Depending  on  the 
end  fixity,  Iq  Is: 


*  0.66 

for  fixed  and  free  (cantilever) 

*  0.63 

for  pin  and  pin  (simple  beam) 

»  0.58 

for  fixed  and  pin 

-  0.52 

for  fixed  and  fixed 

Step  2:  Assume  I/I^  >  1.0  and  calculate  the  amplitude. 
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step  3:  Obtain  a  new  value  of  I/I^  based  on  Initial  ai^illtude. 
Step  4:  Recompute  the  an^lltude  based  on  the  new  value  of  I/I^. 
Step  5:  Repeat  Steps  3  and  4  until  convergence. 

D.5.  STRESSES  DUE  TO  VORTEX  SHEDDING 


Once  the  amplitude  of  vibration  has  been  calculated,  stresses  can  be 
computed  according  to  the  support  conditions.  For  a  simply  supported 
beam  with  a  uniform  load  w,  the  midspan  deflection  Y,  and  the  midspan 
bending  stress  f^^  are  given  as  follows: 


Y 

I 

m 

5  w  1 

15T  *  Tr~  *  ^ 

w 

- 

384  EIY 

%iax 

9 

^bmax 

= 

M  d  .  o  EBY 

T  .  5  -  4.8 

at  midspan 

Expressing 

S 

EOY 

The  K  value  varies 

with 

support  conditions 

and  location  as  shown  on 

Table  0-1. 


Fixity 

Mid-Span 

Ends 

Fix 

Fix 

8.0 

16.0 

Fix 

Pin 

6.5 

11.6 

Pin 

Pin 

4.8 

0 

Fix 

Free 

N.A. 

2.0 

Table  0-1  K  Values  Based  on  Fixity  and  Location 
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The  vortex  shedding  bending  stress  Is  confined  with  the  aeaber  axial 
and  bending  stresses  due  to  global  deformation  of  the  marine 
structure. 

D.6.  FATIGUE  LIFE  EVALUATION 


The  fatigue  life  evaluation  can  be  carried  out  In  a  conservative  two- 
step  process.  First,  the  fatigue  damage  due  to  the  vortex- Induced 
oscillations  Is  calculated  as  0]^.  Second,  a  deterministic  fatigue 
analysis  Is  performed  by  computer  analysis.  Hot  spot  stress  range  vs 
wave  height  (or  wind  velocity)  for  the  loading  directions  considered 
Is  determined  from  the  computer  analysis.  The  critical  direction  Is 
determined  and  a  plot  1s  made.  From  the  plot  of  hot  spot  stress 
range  vs  wave  height  (or  wind  velocity),  the  stress  ranges  for  the 
fatigue  waves  are  determined.  The  maximum  vortex- Induced  stress 
ranges  for  the  fatigue  environment  are  added  to  the  deterministic 
fatigue  stress  ranges.  Then,  the  standard  deterministic  fatigue 
analysis  Is  performed  using  the  Increased  stress  range.  The  fatigue 
damage  calculated  In  this  second  step  is  03.  Therefore  the  total 
fatigue  damage  1s  equal  to  the  sum  of  and  0^,  or  0^^^  «  0^  O^. 

The  fatigue  life  in  years  Is  therefore  calculated  as  I/D^q^. 

A  typical  fatigue  life  evaluation  procedure  Is  given  below: 

Step  1: 

a.  Calculate  the  natural  frequency  fp  (H^)  of  the  member. 

b.  Calculate  the  stability  parameter  of  the  member. 


c.  Determine  the  minimum  required  for  vibrations  based  on  In 
Figure  0-1. 

d.  Calculate  the  minimum  velocity  at  which  current-  or  wind- 
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vortex  shedding  will  occur,  I.e.,  Vm.jp  *  V(req'd)  ^  ^ 

e.  Check  the  applied  velocity  profile  to  see  If  '^max  greater 

than  Vm^p.  If  V^^j^  Is  less  than  V^^p,  then  no  vortex 

oscillations  can  occur. 

f.  For  Vmjix  greater  than  V^^p.  vortex  oscillations  can  occur.  The 

displacement  amplitude  Is  based  on  stability  parameter  K^,  and 
Is  determined  from  Figure  0-2  for  In-line  vibration.  A 

conservative  approach  Is  used  to  determine  Y/d  vs  K^.  For  < 

0.6  the  first  Instability  region  curve  Is  used.  For  >  0.6 
the  second  Instability  region  curve  Is  used.  This 

conservatively  represents  an  envelope  of  maximum  values  of  Y/d 
vs  Kj  from  Figure  0-2.  Olsplacement  amplitude  Is  normalized  to 
Y/d. 

g.  Given  (Y/d),  calculate  the  bending  stress,  f|j. 

h.  Multiply  bending  stress  fjj  by  an  SCF  of  1.5  to  produce  hot  spot 
stress  fjj.  A  larger  SCF  will  be  used  where  necessary. 

I.  From  the  maximum  hot  spot  stress,  the  hot  spot  stress  range  Is 
calculated  as 

J.  Allowable  number  of  cycles  to  failure  (N)  should  be  calculated 
using  an  applicable  S-N  curve  (based  on  weld  type  and 
environment). 

k.  Assume  conditions  conducive  to  resonant  vortex  shedding  occur 
for  a  total  time  of  T  (seconds)  per  annum  (based  on  current  or 
wind  data  relevant  to  applicable  loading  condition). 

l.  Hence,  In  time  T,  number  of  cycles  n  -  fpT  and  the  cumulative 
damage  01  »  n/N  ■  ^pT/N  In  one  year. 

Step  2: 
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a.  Depending  on  marine  structure  In  service  conditions  (I.e. 
structure  In  water  or  In  air)  run  an  applicable  loading 
analysis.  Assuming  a  marine  environment,  run  a  storm  wave 
deterministic  fatigue  analysis  and  obtain  the  results  of  hot 
spot  stress  range  vs  wave  height  for  the  wave  directions 
considered  and  as  many  hot  spots  as  are  needed. 

b.  Determine  the  critical  hot  spot  and  wave  direction  and  draw  the 
hot  spot  stress  range  vs  wave  height  graph. 

c.  Determine  the  hot  spot  stress  range  for  each  of  the  fatigue 
waves. 

d.  For  the  larger  fatigue  waves  In  which  vortex- Induced 
oscillations  occur,  add  the  Increase  In  stress  range  due  to 
vortex- Induced  oscillations  to  the  stress  range  from  the 
deterministic  fatigue  analysis. 

e.  Calculate  the  fatigue  damage  D2  over  a  1  yr  period  for  the  full 
range  of  wave  heights: 

H  n 

“2  •  '  (S) 

f.  Calculate  the  total  fatigue  damage: 

hot  •  Di  +  O2 

g.  Calculate  the  fatigue  life  In  years  as: 
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h.  The  fatigue  life  may  be  modified  to  Include  the  effects  of 
corrosion  pitting  In  caissons.  Corrosion  pitting  produces  an 
SCF  at  the  location  of  the  pit.  The  SCF  Is  calculated  as: 


1  ,  "  W 

SCF  »  - + - r-y 

(1-  |)  (1  - 

where: 

b 

*  pit  depth 

t 

>  nominal  caisson  thickness 

The  new  life  Including  corrosion  damage  Is  calculated  as: 


New  Life  « 


Old  Life 


(SCF) 


7 


This  estimate  of  fatigue  damage  can.  If  necessary,  be  refined  by 
consideration  of  the  nurtjer  of  wave  occurrences  for  different 
directions  and  evaluation  of  the  damage  at  a  number  of  points 
around  the  circumference  of  the  member. 


D.7.  EXAMPLE  PROBLEMS 


0.7.1  Avoidance  of  Wind-Induced  Cross-Flow  Vortex  Shedding 

It  can  be  shown  that  for  a  steel  beam  of  circular  cross  section,  the 
following  relationship  holds: 


where: 


V 


cr 


c. 


Ve 


(L/d) 


7 


'w^  +  wi' 


Vj.y.  -  critical  wind  velocity  of  the  tubular  necessary 
for  the  onset  of  cross-flow  wind-induced  vortex 
shedding 
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c  -  constant  (See  NOTE,  next  page) 

■  reduced  velocity 

Og  >  menber  end  efficiency 

«  I. 5  fixed  ends 

>  1.0  pinned  ends 

Wg  >  weight  per  unit  length  of  tubular 

W]^  >  weight  per  unit  length  of  supported  Itea  (e.g., 

anodes) 

L  «  bean  length 

d  «  tubular  mean  diameter 


For  «  4.7,  Og  »  1.5  (fixed  condition),  and  wj  »  0,  this  reduces 
to: 

»  97240/(L/d)?  ft/sec 
»  29610/(L/d)'^  m/sec 

Hence,  If  maximum  expected  wind  speed  Is  65.6  ft/s  (20  m/s),  then 
setting  all  brace  L/d  ratios  at  38  or  less  precludes  wind-induced 
cross-flow  vortex  shedding,  and  no  further  analyses  or  precautions 
are  required. 

However,  maximum  wind  speeds  may  be  so  high  that  the  above  approach 
may  be  uneconomical.  In  this  case,  either  precautionary  measures 
must  be  taken  or  additional  analyses  considering  strength  and  fatigue 
must  be  undertaken. 
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NOTE: 

The  relationship  given  Is  based  on: 

a  4  % 

''cr  ■  V  'n  '*  ■  V  (  5T  •  " 

substituting 

»n  “  <"e 

mi  »  (Wq  +  w^)/g 
I  -  ir  d^t/8 

E  -  4176  X  10®  Ibs/ft^  (200,000  MH/m^) 
g  -  32.2  ft/sec^  (9.806  m/sec^) 


Wq  -  Yj  «dt 

Yj  *  weight  density  of  steel,  490  1bs/ft^  (0.077  MN/ra^) 
2  2 


''cr’  V  (“V— )  I  I  l’*  •  d 

(''o  ^  ^'l) 

E  (w  /y.)  g  . 

V  V  (-r-)  I  — — — r 


Substituting  for  E,  y-.  and  g 


cr 

where 


(L/d) 


'«  I  " 

ITT  I  ^ 


0  '^1 


constant  C  *  9195 
-  2800 


for  as  ft/sec 

for  as  m/sec 
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Analysis  for  Wind-Induced  Cross-flow  Vortex  Shedding 

Using  procedures  discussed  in  Section  0.4  a  flare  structure  bracing 
meoibers  are  ana1y2ed  for  crossflow  oscillatons  produced  by  vortex 
shedding.  The  analysis  is  performed  using  a  Lotus  spreadsheet.  The 
general  procedure  is  as  loliows: 

(a)  Member  and  environmental  parameters  are  input. 

(b)  Critical  velocity,  peak  amplitudes  of  oscillation  and 
corresponding  stress  amplitudes  are  coii4)uted. 

(c)  The  time  (in  hours)  of  crossflow  oscillation  required  to  cause 
fatigue  failure  is  computed. 

Analysis  Description 

The  following  is  a  detailed  description  of  the  spread  sheet  input  and 
calculation. 

(a)  Spread  Sheet  Terminology 

Columns  are  labeled  alphabetically  while  rows  are  labeled 
numerically.  A  "cell"  is  identified  by  referring  to  a  specific 
row  and  column. 

(b)  General  Parameters 


The  following  are  parameters  common  to  all  members  analyzed  as 
given  at  the  top  of  the  spread  sheet. 

CELL  C5:  DAMPING  RATIO  -  e 

CELL  C6:  AIR  MASS  DENSITY  »  p 
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CELL  C7:  KINEMATIC  VISCOSITY  -  u 

CELL  C8:  STRESS  CONCENTRATION  FACTOR  »  SCF 

CELL  C9:  MODULUS  OF  ELASTICITY  -  E 

CELL  K5:  RATIO  OF  GENERALIZED  MASS  TO  EFFECTIVE  MASS  «  (  ) 

e 

CELL  K6:  FIXITY  PARAMETER  IN  FORMULA  FOR  CRITICAL  VELOCITY  -  n 

e 

CELL  K7:  FIXITY  PARAMETER  IN  FORMULA  FOR  STRESS  AMPLITUDE  «  C 
CELL  K8:  FIXITY  PARAMETER  IN  FORMULA  FOR  MEMBER  FREQUENCY  > 

(c)  Specific  Member  Analysis 

The  following  describes  the  content  of  each  column  In  analyzing 
a  specific  member.  Entries  and  formulas  for  vortex  shedding 
analysis  of  member  group  HI  on  line  16  are  also  provided. 
Formula  coding  Is  described  In  the  LOTUS  1-2-3  Users  Manual. 

COLUMN  A:  ENTER  THE  MEMBER  GROUP  IDENTIFIER 

COLUMN  B:  ENTER  THE  EFFECTIVE  SPAN  OF  THE  MEMBER  -  L  (m) 

COLUMN  C:  ENTER  THE  OUTSIDE  DIAMETER  OF  THE  TUBULAR  »  d  (mm) 

COLUMN  D:  ENTER  THE  TOTAL  OUTSIDE  DIAMETER  «  D  (mm)  INCLUDING  AS 
APPLICABLE.  MARINE  GROWTH.  FIRE  PROTECTION.  ETC. 

COLUMN  E:  ENTER  THE  TUBULAR  WALL  THICKNESS  -  t  (mm) 

COLUMN  F;  ENTER  ADDED  MASS  (kg/m).  IF  APPLICABLE 

COLUMN  6:  THE  MOMENT  OF  INERTIA  OF  THE  TUBULAR  -  I  (cm^)  IS 
COMPUTED. 


D-21 


I  -  -J- 1  (-f  - 


COLUMN  H:  THE  TOTAL  EFFECTIVE  MASS  IS  COMPUTED 


"•e  *  “r  ■  "a 


COLUMN  I:  THE  CRITICAL  VELOCITY  FOR  CROSSFLOW  OSCILLATION  IS 
COMPUTED. 


cr 


13160  n‘ 
(L/0)2 


(«/s) 


COLUMN  J:  ENTER  THE  THRESHOLD  WIND  VELOCITY  -  V^^,^ 


COLUMN  K:  THE  STABILITY  PARAMETER  IS  COMPUTED 

2«.  (2.E) 

K  «  — ^ - 

*  PO^ 

COLUMN  L:  THE  REYNOLDS  NUMBER  IS  COMPUTED 


Before  performing  calculation  In  the  following  columns,  the  critical 
velocity  Is  con^ared  with  the  threshold  velocity.  If  the  critical 
velocity  Is  larger,  crossflow  oscillations  will  not  occur  and  the 
computations  are  supressed.  An  "N.A."  Is  then  Inserted  In  each 
column. 


If  the  critical  velocity  Is  less  than  the  threshold  value,  the 
following  computations  are  performed. 


COLUMN  M:  THE  AMPLITUDE  OF  VIBRATION  IS  COMPUTED 


V  _  g  D  Where  g  ■  0.04925  for  Re  >  500,000 

/mJ  %  ^  and  g  ■  0.0717B  for  Re  <  500,000 

'me  ’  s 
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COLUMN  N:  THE  STRESS  AMPLITUDE  IS  COMPUTED 


fb  -  (MPa) 

where  C  depends  on  beam  end  fixity  (see  Section  D.4) 

COLUMN  0:  THE  HOT  SPOT  STRESS  RANGE  IS  COMPUTED 

S  -  2  (SCF)  fjj  (MPa) 

COLUMN  P:  THE  NUMBER  OF  CYCLES  TO  FAILURE  UNDER  THE  HOT  SPOT 
STRESS  RANGE  IS  COMPUTED. 

H  •  ‘•‘>910^)  (cycles) 


COLUMN  Q:  THE  MEMBER  NATURAL  FREQUENCY  IS  COMPUTED 


'n- 


*n  ,  El 

rr  ( —3 


(Hj) 


where  a^  depends  on  beam  end  fixity  (see  Section  0.2) 


COLUMN  R:  THE  TIME  IN  HOURS  TO  FATIGUE  FAILURE  UNDER  N  CYCLES  OF 
STRESS  RANGE  S  IS  COMPUTED 


0.8.  METHODS  OF  MINIMIZING  VORTEX  SHEDDING  OSCILUTIONS 


D.8.1  Control  of  Structural  Design 

The  properties  of  the  structure  can  be  chosen  to  ensure  that  critical 
velocity  values  In  steady  flow  do  not  produce  detrimental 
oscillations. 
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Experiments  have  shown  that  for  a  constant  mass 

2 

parameter  (m/pd  =  2.0),  the  critical  velocity  depends  mainly  on  the 
submerged  length/diameter  (L/d)  ratio  of  the  member. 

Thus,  either  high  natural  frequency  or  large  diameter  Is  required  to 
avoid  VIV's  In  quickly  flowing  fluid.  A  higher  frequency  will  be 
obtained  by  using  larger  diameter  tubes,  so  a  double  benefit 
occurs.  An  alternative  method  of  Increasing  the  frequency  Is  to 
brace  the  structure  with  guy  wires. 

D.8.2  Mass  and  Damping 

2 

Increasing  the  mass  parameter,  m/pd  ,  and/or  the  damping  parameter 
reduces  the  amplitude  of  oscillations;  If  the  Increase  Is  large 
enough,  the  motion  Is  suppressed  completely.  While  high  mass  and 
damping  are  the  factors  that  prevent  most  existing  structures  from 
vibrating,  no  suitable  design  criteria  are  presently  available  for 
these  factors,  and  their  effects  have  not  been  studied  In  detail. 

Increasing  the  mass  of  a  structure  to  reduce  oscillatory  effects  may 
not  be  entirely  beneficial.  The  Increase  may  produce  a  reduction  In 
the  natural  frequency  (and  hence  the  flow  speeds  at  which  oscillation 
will  tend  to  occur).  It  Is  thus  possible  that  the  addition  of  mass 
may  reduce  the  critical  speed  to  within  the  actual  speed  range. 
However,  If  increased  mass  is  chosen  as  a  method  of  limiting  the 
amplitude  of  oscillation,  this  mass  should  be  under  stress  during  the 
motion.  If  so,  the  mass  will  also  contribute  to  the  structural 
damping.  An  unstressed  mass  will  not  be  so  effective. 

If  the  structure  Is  almost  at  the  critical  value  of  the  combined 
mass/damping  parameter  for  the  suppression  of  motion,  then  a  small 
additional  amount  of  damping  may  be  sufficient. 
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D.8.3  Devices  and  Spoilers 

Devices  that  modify  flow  and  reduce  excitation  can  be  fitted  to 
tubular  structures.  These  devices  (see  Figure  D-5)  work  well  for 
isolated  members  but  are  less  effective  for  an  array  of  piles  or 
cylinders.  Unfortunately,  there  1s  no  relevant  Information 
describing  how  the  governing  stability  criteria  are  modified.  The 
most  widely  used  devices  are  described  below. 

Guy  Wires 

Appropriately  placed  guy  wires  may  be  used  to  Increase  member 
stiffness  and  preclude  wind-induced  oscillations.  Guy  wires  should 
be  of  sufficient  number  and  direction  to  adequately  brace  the  tubular 
member;  otherwise,  oscillations  may  not  be  eliminated  completely  and 
additional  oscillations  of  the  guys  themselves  may  occur. 

Strakes  or  Spoilers 

Strakes  and  spoilers  consist  of  a  number  (usually  three)  of  fins 
wound  as  a  helix  around  the  tubular.  These  have  proven  effective  In 
preventing  wind-induced  cross-flow  oscillations  of  structures,  and 
there  is  no  reason  to  doubt  their  ability  to  suppress  in-line  motion, 
provided  that  the  optimum  stroke  design  is  used.  This  comprises  a 
three-star  helix,  having  a  pitch  equal  to  five  times  the  member 
diameter.  Typically  each  helix  protrudes  one-tenth  of  the  member 
diameter  from  the  cylinder  surface.  To  prevent  In-line  motion, 
strakes  need  only  be  applied  over  approximately  In  the  middle  one- 
third  of  the  length  of  the  tubular  with  the  greatest  amplitude. 
Elimination  of  the  much  more  violent  cross-flow  motion  requires  a 
longer  strake,  perhaps  covering  the  complete  length  of  tube.  The 
main  disadvantage  of  strakes,  apart  from  construction  difficulties 
and  problems  associated  with  erosion  or  marine  growth,  is  that  they 
increase  the  time-averaged  drag  force  produced  by  the  flow.  The  drag 
coefficient  of  the  straked  part  of  the  tube  is  independent  of  the 
Reynold's  number  and  has  a  value  of  Cq  «  1.3  based  on  the  tubular 
diameter. 


D-25 


Shrouds 


Shrouds  consist  of  an  outer  shell,  separated  from  the  tubular  by  a 
gap  of  about  0.10  diameter,  with  many  small  rectangular  holes.  The 
limited  data  available  Indicates  that  shrouds  may  not  always  be 
effective.  The  advantage  of  shrouds  over  strakes  Is  that  their  drag 
penalty  Is  not  as  great;  for  all  Reynold's  numbers,  Cq  «  0.9  based  on 
the  Inner  tubular  diameter.  Like  strakes,  shrouds  can  eliminate  the 
In-line  motion  of  the  two  low-speed  peaks  without  covering  the 
complete  length  of  the  tubular.  However,  any  design  that  requires 
shrouds  (or  strakes)  to  prevent  cross-flow  motion  should  be 
considered  with  great  caution.  Their  effectiveness  can  be  minimized 
by  marine  growth. 

Offset  Dorsal  Fins 


This  Is  the  simplest  device  for  the  prevention  of  oscillations.  It 
Is  probably  the  only  device  that  can  be  relied  upon  to  continue  to 
work  In  the  marine  environment  over  a  long  period  of  time  without 
being  affected  adversely  by  marine  growths  It  has  some  drag  penalty, 
but  this  Is  not  likely  to  be  significant  for  most  designs. 

The  offset  dorsal  fin  Is  limited  to  tubular  structures  that  are 
subject  to  In-line  motion  due  to  flow  from  one  direction  only  (or  one 
direction  and  its  reversal,  as  In  tidal  flow). 

This  patented  device  comprises  a  small  fin  running  down  the  length  of 
the  tubular.  Along  with  the  small  drag  Increase  there  1s  a  steady 
side  force.  This  may  be  eliminated  In  the  case  of  the  total  force  on 
multi-tubular  design  by  placing  the  f1n  alternately  on  opposite  sides 
of  the  tubulars. 
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Figure  0'*4  Asplitude  of  Response  for  Cross-Flow  Vibrations 
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Figure  D-5  Typical  Devices  and  Spoilers 
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